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ABSTRACT 

We examine the temperature history of gas accreted by forming galaxies in smoothed 
particle hydrodynamics (SPH) simulations. About half of the gas follows the track 
expected in the conventional picture of galaxy formation, shock heating to roughly 
the virial temperature of the galaxy potential well (T ^ 10^ K for a Milky Way type 
galaxy) before cooling, condensing, and forming stars. However, the other half radiates 

■ its acquired gravitational energy at much lower temperatures, typically T < 10^ K, and 

the histogram of maximum gas temperatures is clearly bimodal. The "cold mode" of 

, gas accretion dominates for low mass galaxies (baryonic mass Mgai < 10^'^-^ Mq or halo 

' mass Mhaio ^ 10^^ '* M©), while the conventional "hot mode" dominates the growth of 

high mass systems. Cold accretion is often directed along filaments, allowing galaxies 
to efficiently draw gas from large distances, while hot accretion is quasi-spherical. The 
galaxy and halo mass dependence leads to redshift and environment dependence of 
• cold and hot accretion rates, with cold mode dominating at high redshift and in low 

[ density regions today, and hot mode dominating in group and cluster environments at 

. low redshift. The simulations reproduce an important feature of the observed relation 

' between galaxy star formation rate (SFR) and environment, namely a break in star 

Q^, formation rates at surface densities E ~ Ift-fgMpc"^, outside the virial radii of large 

I ' groups and clusters. The cosmic SFR tracks the overall history of gas accretion, and 

O ' its decline at low redshift follows the combined decline of cold and hot accretion rates. 

The drop in cold accretion is driven by the decreasing infall rate onto halos, while for 
^ ' hot accretion this slower mass growth is further modified by the longer cooling times 

' within halos. If we allowed hot accretion to be suppressed by conduction or AGN 

feedback, then the simulation predictions would change in interesting ways, perhaps 
resolving conflicts with the colours of ellipticals and the cutoff of the galaxy luminosity 
^ . , function. The transition at A/haio ~ 10^^ *Mq between cold mode domination and 

■ hot mode domination is similar to that found bv lBirnboim fc Dekell pOOS.'l using 1-d 

simulations and analytic arg uments. The corre s pondi ng baryonic mass is tantalisingly 
close to the scale at which iKauffmann et all l)2003l) find a marked shift in galaxy 
properties, and we speculate on possible connections between these theoretical and 
observational transitions. 
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1 INTRODUCTION 

The conventional sketch of galaxy formation has its roots 
in classic papers of the late '70s and early '80s, with 
initial discussio n s of colla pse and co oling criteria by 
iRees fc Ostriked Jl977|) and ISilkl (ll977D, t he addition of 
dark matter halos bv Iwhite fcRee^ Jl978l). a nd the disk 
formation model of iFal^m^Efete^iiounilQSOh . According 
to this sketch, gas falling into a dark matter potential 
well is shock heated to approximately the halo virial tem- 



perature, Tvir = 10^(t;circ/167kms~^)^ K, putting it in 
quasi-hydrostatic equilibrium with the dark matter. Gas 
in the dense, inner regions of this shock heated halo ra- 
diates its thermal energy, loses its pressure support, set- 
tles into a centrifugally supported disk, and forms stars. 
Mergers of disks can scatter stars onto disordered orbits, 
producing spheroidal systems, which may regrow disks if 
they experience subsequent gas accretion. Over the last 
decade, the ideas of these seminal papers have been up- 
dated and extended into a powerful "semi-analytic" frame- 
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work for galaxy formation calculations (e. g., 'White fc Frenk* 
199l|: Kauffmann. White^ fc Guiderdonil |1993: Cole et al., 
I994I: Uvila-Reese et alJ Il99d IMo. Mao, fc Whitd Il998t 
Somerville fc Primacklll999t) . 

The geometry seen in N-body and hydrodynamic cos- 
mological simulations, where the densest structures often 
form at the nodes of a filamentary network, is clearly more 
complicated than the spherical geometry underlying semi- 
analytic calculations. Nonetheless, a substantial fraction of 
the gas in these simulations does shock heat to T ~ Tvir, 
and some of this gas does cool and settle into galaxies. 
The approximate agreement between semi-analytic models 
and smoothed parti cle hydrodynamics (SPH) calculations of 
galaxy m asses (e.g.. 'Benson et alJl200ll:lYoshida et alJl2002l: 
iHelly et ~al. 2003) has therefore been taken as evidence that 
the conventional sketch, while idealised, captures most of 
the essential physics. In this paper, we use SPH simulations 
of cosmological volumes to argue that this sketch requires 
an important revision: roughly half of the gas accreted by 
the simulated galaxies is never shock heated close to the 
halo virial temperature (T ~ 10^ K for a Milky Way type 
galaxy) but instead radiates its acquired gravitational en- 
ergy from r<2.5 X 10^ K (often T<5 x lO^K). This "cold 
mode" of gas accretion dominates for lower mass galaxies 
(baryon mass Mgai < 2 x W^^Mq), while the conventional, 
"hot mode" of gas accretion dominates the growth of high 
mass systems. As a result, "cold mode" accretion dominates 
at high redshift {z > 3) and in low density environments to- 
day, while "hot mode" accretion dominates in group and 
cluster environments at low redshift. 

There is, in fact, a long history of results suggesting 
that cold accretion could be an important element of galaxy 
formation. Binney (1977), using analytic models of proto- 
galaxy collapse, argued that the amount of shock heating 
could be small for plausible physical conditions, with only 
a fraction of the gas reaching temperatures T ~ Tvir. In 
the first SPH simulations of forming galaxies (Katz & Gunn 
1991), which had idealised initial conditions but included 
small scale power leading to hierarchical formation, most 
of the gas never heated above T ~ 3 x 10^ K, with much 
of the cooling radiati on t herefore emerging in the Lya line. 
iKatz fc White (1993^) and' Katz et alj il994h showed the im- 
portance of filamentary structures as channels for gas ac- 
cretion in simulations with cold dark matter (CDM) initial 
conditions. Recent studies based on SPH simulations of cos- 
mological volumes reveal the situation even more starkly. 
Fardal et al. (2001) showed that most of the cooling radia- 
tion in their simulations comes from gas with T < 2 x 10^ K, 
again implying that a significant fraction emerges in the Lya 
line. Since gas starting at T ~ 10® K must radiate 90% 
of its therm a l ener gy by the time it cools to T ~ 10^ K, 
iFardal et al.l (I2OOII) concluded that the majority of the gas 
entering galaxies (indeed, the majority of the gas experienc- 
ing any significant cooling) must not be heated to the virial 
temperature of any dark matter halo resolved by the sim- 
ulation. Kay et al. (2000) directly tracked the temperature 
histories of particles that ended up in their simulated galax- 
ies and found that only 11% of these particles were ever 
heated to a temperature above 10^ K. 

Motivated by these results, we here investigate the tem- 
perature histories of accreted gas particles as a function of 
galaxy mass, redshift, and environment, thoroughly quanti- 



fying the relative importance of the cold and hot modes of 
gas accretion in our SPH simulations. We use several simu- 
lations to demonstrate the insensitivity of our primary con- 
clusions to numerical resolution over a wide dyn amic range. 
We re p orted initial result s from our study in iKatz et alJ 
J2003h . Ifiirnboim fc Dekell ||200 J) have investigated similar 
issues with a complementary approach based on high resolu- 
tion, spherically symmetric collapse calculations. They find 
that a virial shock fails to develop if the gas cooling time 
is shorter than the local dynamical time, so in these cases 
gas shells fall far inside the halo virial radius without ever 
being heated to high temperature. iBirnboim fc DekeJll2003h 
show that the cooling time condition corresponds approxi- 
mately to a threshold in the galaxy's halo mass, with little 
dependence on redshift. As we will show in fj^below, our re- 
sults are in quite good agreement with Birnboim fc Dekel's, 
despite the radically different approaches. Traditional semi- 
analytic models also distinguish between halos with rapid 
post-sh ock co oling and halos with slow post-shock cooling 
(e.g.. iWhite fc Frcnk 1991') . though the distinction has re- 
ceived relatively little attention in discussions of these mod- 
els, and heating to the virial temperature is assumed in ei- 
ther case. In quantitative terms, we find that cold accre- 
tion plays a much larger role in our simulations than rapid- 
cooling accretion plays in standard semi-analytic calcula- 
tions (see i|6.3l and the Appendix). 

In addition to being an important aspect of the physics 
of galaxy formation, the existence of distinct cold and 
hot modes of gas a ccretion could have interesting observa- 
tional implications. iFardal et al.l (1200111 emphasised one of 
these implications: cold accretion allows much of the cool- 
ing radiation associated with galaxy formation to emerge 
in the Lyg Une instead of the X- ray continuum (see also 
iHaiman. Soaans. fc Quataertll2o'ool) . A reduced role for hot 
accretion might help explain why diffuse X-ray emission 
from late-type galaxy halos is well b elow the prediction s of 
standard semi-analytic calculations jBenson et al.ll2bo O'). A 
second class of implications relat es to the cosmic star forma- 
tion history. Murali et al. J2002t hereafter MKHWD) show 
that galaxies in SPH simulations, like the ones analysed here, 
gain most of their mass through smooth accretion of gas, 
not through mergers with pre-existing galaxies (at least not 
galaxies above the simulation's resolution limit). They fur- 
ther find that the global history of star formation tracks 
the global history of gas accretion rather than the merger 
history. In subsequent analysis, we have found that these 
generalisations hold fairly well on a galaxy-by-galaxy basis 
(Mailer et al., in preparation). Thus, understanding gas ac- 
cretion is nearly tantamount to understanding the history 
of star formation, at least in the simulations. 

The cosmic star formation history inferred from near- 
UV luminos ity functions decline s sharply between z ~ 1 and 
2 = (e.g.. lMadau et alJll99d) . While semi-analytic mod- 
els and hydrodynamic simulations both predict a drop in 
the global star formation rate (SFR) over this redshift in- 
terval, it is difficult to explain the full order-of-ma gnitude 
reduction implied by the data jBaugh et al.l 12004 . With- 
out an ad hoc fix, semi-analytic and numerical calculations 
also predict continuing gas accretion and star formation in 
old, massive galaxies, so they have difficulty reproducing the 
"red envelope" of th e observed galaxy population as a func- 
tion of redshift (e.g.. lCole et aLlboOO) . Of course, partition- 
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ing the accretion into cold and hot modes does not, in itself, 
change the simulation predictions, but it may illuminate the 
physics behind the predicted drop in the SFR and explain 
differences between numerical and semi-analytic results. The 
declining SFR is frequently attributed to the longer cool- 
ing times in the hotter, lower de nsity halos that prevail at 
low redshift fe.g.. iBlanton et aljf200fl ). but the efficiency of 
cold accretion could be strongly affected by other factors. 
Furthermore, stellar or AGN feedback, or heat conduction, 
could have different effects on the cold and hot accretion 
modes, since the incoming gas has different geometry and 
density. Allowing for such a difference in theoretical models 
could have interesting observational consequences, since the 
relative importance of cold and hot accretion depends on 
redshift, environment, and galaxy mass. 



The sensitivity of cold and hot accretion rates to envi- 
ronment could also play a role in e xplaining the well known 
morphology-density relation (e.g., iHubbld Il936l : iDresslerH 
Ll980; Postman fc G eller f 98i) and the associated correla- 
tion between galaxy SFR and local density (e.g.,[ 



et alJ 



120021: iGomez et aLll2003l : iKauffmann et alj|2004) . Mergers, 
ram pressure stripping, ll^^^^^Go^^72l ). galaxy "ha- 
rassment" by weak perturbations in clusters (iMoore et alJ 
[l996; Moore. Lake, fc Katz. 1998.). truncation of gas supplies 
Tinslev fc Caldwell I 1198(1: ISomerviUe fc Primackl 
[ 52^)!^and longer c ooling times in hotter environments 
Irnlanton et alJ I2OOOI) may all contribute to the origin of 
these correlations. However, environmental effects that shut 
off cold accretion flows could suppress star formation and 
disk growth, and they might explain why transitions in 
galaxy properties appear t o start well beyond t he virial 
radii o f groups and cluste r s (iLewis et alj E002: Gom ez et all 
l2003l) . IKauffmann et alJ ||2003) find a clear transition in 
galaxy properties at a baryonic mass Mb ~ 3 x IO^^A/q, 
with lower mass galaxies having active star formation, low 
surface mass density, and a disk morphology while higher 
mass galaxies have old stellar populations, high surface mass 
density, and a bulge-dominated morphology. The transition 
that we find between galaxies dominated by cold accretion 
and galaxies dominated by hot accretion occurs at a similar 
mass scale and could be connected to this broader transition 
in galaxy properties. 



After briefly describing our simulations and analysis 
methods in ^ we present our basic results on the global 
signiflcance and mass dependence of cold and hot accretion 
in fJ3| We assess the numerical robustness of these results 
in in particular comparing different simulations to show 
that cold mode gas accretion is found in simulations with a 
large dynamic range in mass resolution. We investigate the 
dependence of the cold and hot accretion rates and the corre- 
sponding star formation rates on galaxy environment in ^ 
including a comparison of the predicted correlation s between 
SFR and environment to the lGomez et al] ll2003f) observa- 
tions. In f|S|we describe our current physical understanding 
of the cold and hot accretion modes and discuss some of the 
potential implications mentioned above. We summarise our 
conclusions in iJ7| 



2 SIMULATIONS AND NUMERICAL 
METHODS 

2.1 Simulation parameters 

We adopt an inflationary cold dark matter model domi- 
nated by a cosmological constant, ACDM, with Q.m = 0.4, 

= 0.6, h = Ho/{im km s"^ Mpc"^) = 0.65, and 
a primordial power spectrum index n — 0.93. For the 
amplitude of mass fluctuations we use erg = 0.8, which for 
our adopted parameters is consisten t both with COBE 
norm ali sation using CMBFAST dSeliak fc Zaldariaeal 
il996l: IZaldarriaea. Seliak fc Bertschinge j ll998D and 
with the o bserved abundance of rich clusters 
iWhite. Efstath iou. fc Fronk 1993). For the baryonic 
density we adopt fii, — 0.02/i~^, a value consistent both 
with the deu terium abundance in high redshift Lyman 
limit systems jBurles fc Tvtlerlll998h and the value derived 
from cosmic microwave ba ckground (CMB) anisotropy 
measurements llde Bernardis c t al. 200J.). Our values of 
Qb, n, and as are close to those inferred by recent joint 
analysis of CMB an isotropy measurements from WMAP 
iBennett et al J 120031) and galax;y clust ering data from the 
2dF Galaxy Redshift Survey (2dFGRS; Colloss ot al. 20^ 
and the Sloan Digital Sky Survey (SDSS;,York ct al. 200j3). 
while our assumed f2m is higher by about 1.5a. We have 
recently repeated one of our runs using t he parameter 
values implied by the WMAP analysis of ISoergel et al.l 
i2003l) . and preliminary investigation shows results similar 
to those reported here. 

Our primary results are derived from a simulation that 
models a 22.222/i~^ Mpc comoving periodic cube using 128'^ 
dark matter particles and 128'' gas particles. Gravitational 
forces are softened using a cubic spline kernel of comoving 
radius 5h~^ kpc, approximately equivalent to a Plummer 



force softening of eg 



3.5/i kpc. Our baryonic mass 



threshold for resolved galaxies (see is 6.8 x W^Mq, the 
mass of 64 gas particles, and there are 1120 galaxies in the 
box above this threshold ai z = 0. To approximately match 
our galaxy mass scale to an obser ved luminosity scale, we 
note that the Bla nton et all i2003h r-band luminosity func- 
tion yields a space density of 0.0032 Mpc~^ for galaxies 
bright er than the characteristic luminosity L* of a iSchechteil 
||l97g) function fit. The baryonic mass threshold that yields 
the same space density in our simulation is 2.45 x W^^ Mq, 
so if we identify this mass with L« and assume that lumi- 
nosity is approximately proportional to mass, our resolution 
threshold corresponds roughly to L»/36. 

We draw on six additional simulations to investigate the 
influence of mass resolution and the presence of a UV back- 
ground fleld on our results. Parameters of all the simulations 
are listed in Table Q A simulation of comoving box length 
X Mpc and A*'^ particles is designated Lx/N; thus, our 
primary simulation is L22/128. Most of the simulations were 
run with a photoionizing UV background (see below); those 
that were not are designated "nb" for "no background." The 
fuU suite of simulations — L50/144nb, L22/128, L22/64nb, 
Lll/64, Lll/64nb, Lll/128, L5. 5/128 — spans a factor of 
512 in mass resolution. However, the highest resolution sim- 
ulations (Lll/128 and L5. 5/128) have been evolved only to 
z = 3, and the lowest resolution simulations have no UV 
background. 
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Name 


L(/i-i Mpc) 


N 


2fln 


uv 


Mrcs{MQ) 


L50/144nb 


50 


2 X 144^ 





No 


5.4 X 10^" 


L22/128 


22.22 


2 X 128^ 





Yes 


6.8 X 10'' 


L22/64nb 


22.22 


2 X 643 





No 


5.4 X 10^° 


Lll/64 


11.11 


2 X 64^ 


3 


Yes 


6.8 X lO'* 


Lll/64nb 


11.11 


2 X 64^ 


3 


No 


6.8 X lO'* 


Lll/128 


11.11 


2 X 1283 


3 


Yes 


8.5 X 10** 


L5. 5/128 


5.55 


2 X 128^ 


3 


Yes 


1.1 X 10** 



Table 1. Parameters of the simulations used in this paper. L is the comoving box size, N is the total number of particles (dark+baryonic), 
2fin is the final redshift to which the simulation has been evolved, UV indicates whether or not a UV background is included in the 
calculation of cooling and heating rates, and Afros is the baryonic mass resolution threshold, corresponding to the mass of 64 gas particles. 



2.2 The Simulation Code 

Our simulations are perfor med using the paral- 
lel version of TreeSPH (iHernauist fc Katd Il989l: 
Katz. We inberg, fc Hernauist' 'l99d, hereafter KWH; 
Dave. Du biriski. & Hernauist 1997). This co de combines 
smoothed particle hydro dynamics (SPH; iLucvl 119771: 
iGingold fc MohagharJll997^ with the hierarc hical tree algo- 
rithm for computation of gravitational forces iBarnes fc Hud 
ll986l:lHerna'uislll987^ ■ TreeSPH is a completely Lagrangian 
code, adaptive both in space and in time. In our simu- 
lations, gas properties are estimated by smoothing over 
32 nearby particles. There are three kinds of particles in 
our simulations: dark matter, stars and gas. Colhsionless 
particles (dark matter and stars) are influenced only by 
gravity, while gas particles are influenced by pressure 
gradients and shocks in addition to gravity. 

We use the geometric averaging form of the energy equa- 
tion (Hernauist & Katz 1989). Gas particles experience adi- 
abatic heating and cooling, shock heating, inverse Comp- 
ton cooling off the microwave background, and radiative 
cooling via free-free emission, colhsional ionization and re- 
combination, and collisionally excited line cooling. We as- 
sume primordial abundances (since we are primarily inter- 
ested in following gas before it gets into galaxies), and we 
include only atomic cooling processes, so gas cannot cool 
below T ~ 10" K. In most of our simulations (see above), 
we include photoionization by a spatially uniform UV back- 
ground, which heats low temperature gas and suppresses 
coohng processes involving neutral atoms at low gas densi- 
ties. For the spectral shape and intensity of the UV back- 
ground we use the calculations of Haardt & Madau (1996). 
Our calculations of cooling and heating rates are discussed 
in detail by KWH, who also illustrate the influence of pho- 
toionization on these rates. 

We heuristically include star formation and its associ- 
ated supernova feedback, as described by KWH. In brief, gas 
with physical density pgas > O.lmifcm"^ (hydrogen number 
density uh > O.lcm"^) is assumed to convert into stars on 
a timescale set by the dynamical time or the cooling time, 
whichever is longer. We also require star-forming gas to be 
Jeans unstable, part of a converging flow (V • v < 0), and 
above the virial overdensity (pgas/pbar > 55.7), but gas that 
satisfies the physical density criterion usually satisfies the 
other three criteria as well. Gas reaches this high density 
only after cooling to T ~ 10'* K, and we implicitly include 
the subsequent molecular and metal-line cooling to lower 
temperatures as part of the star formation process. Our for- 



mula for the star formation rate leads to a relation with 
gas surface dens ity similar to a Schmidt law llSchmiddll959l : 
lKennicut3ll998^ . At a technical level, we convert gas parti- 
cles to colhsionless star particles by means of intermediate 
particles that feel reduced gas forces (see KWH). This allows 
us to trace each star particle back to a unique gas progenitor. 

During the formation of stars, supernova feedback en- 
ergy is added to the surrounding gas particles in the form 
of heat. This energy is added gradually with an exponen- 
tial time decay of 2 x 10^ years. We calculate this energy 
assuming that stars with masses above 8Mq explode as su- 
pernovae, which for our adopted Miller-Scalo initial mass 
function (Miller & Scalo 1979) gives 7.35 x 10~^ supernovae 
per solar mass of formed stars. Each supernova deposits 10^* 
ergs of energy into the surrounding medium. The surround- 
ing medium is usually dense, so the deposited energy is typi- 
cally radiated away before it can drive a galactic scale wind. 

2.3 Identifying Galaxies and Hales 

Cosmological simulations that incorporate cooling 
and star formation produce dense groups of bary- 
onic particles wit h the sizes and masses o f ob- 
served galaxies (iKatz. iernauist. fc Weinberg! Il992l : 
lEvrard. Summers, fc Davia ^9^1 ~To identify these dense 
groups we use the group finding alg orithm Spline Kernel 
Inter polative DENMAX (SKID)^ (iGelb fc Bertschingeil 
I1994I : KWH). This algorithm involves four basic steps: 
(1) determine the smoothed baryonic density field; (2) 
move baryonic particles towards higher density along the 
initial gradient of the baryonic density field; (3) define the 
initial group to be the set of particles that aggregate at a 
particular density peak; (4) link together initial groups that 
are very close together; (5) remove particles from the group 
that do not satisfy a negative energy binding criterion 
relative to the group's centre of mass. We apply SKID to 
the population of all star particles and those gas particles 
that have temperatures T < 3 x 10* K and overdensities 
Pgas/pgas > 10'^, and we henceforth refer to the groups of 
stars and cold gas that SKID identifies simply as "galaxies." 
Tests on simulations with varying mass resolution show 
that the simulated galaxy population becomes substan- 
tially incomplete below a baryonic mass corresponding to 
~ 64mspH but is fairly robust above this limit (see, e.g., 
MKHWD). We therefore adopt 64mspH (6.8 x 10^ Mq in the 
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L22/128 simulation) as our resolution threshold and ignore 
lower mass galaxies in our analysis. Because of our high 
overdensity threshold for star formation, essentially all star 
formation in the simulation takes place in galaxies, though 
some of these are below the resolution limit, and some 
stars are tidally stripped from galaxies during dynamical 
interactions. 

We identify d ark matter halos u sing a friends-of-friends 
(FOF) algorithm dOavis et alJll985l) . which selects groups 
of particles in which every particle has at least one neigh- 
bour within a specified linking length. We choose the link- 
ing length to correspond to the interparticle separation at 
1/3 of the virial overdensity pvir/p, which is calculated for 
the value of Q,m at e ach redshift using the fitting formula 
of Rcitavama fc Suto ] 11996). We refine the halos and as- 
sign a virial mass, virial radius, and virial temperature us- 
ing a spherical overdensity (SO) criterion. Specifically, we 
set the centre of the group to be at the position of the 
most bound FOF particle, and we go out in radius until 
the mean enclosed overdensity (dark matter plus baryons) 
equals pvir{z)/p. We define this radius to be 7?vir, the mass 
within it to be Mvir, and the halo circular velocity to be 
^^circ = (GMvir/flvir)^''^. We define the halo virial tempera- 
ture by fcTvir = so it represents the temperature 
at which the gas would be in hydrostatic equilibrium if the 
potential well were isothermal. Since gas at this tempera- 
ture would be fully ionized, we adopt /i = 0.59, appropriate 
to fully ionized, primordial composition gas. We associate 
galaxies with these, refined, SO halos. 

Our minimum galaxy mass Mgai,min = 64mspH cor- 
responds to a minimum host halo mass Mhaio.min ~ 
(flm/f2i))A/gai,min, siucc the fractiou of cold gas in a halo 
never exceeds the universal baryon fraction by a large fac- 
tor (though it can sometimes be slightly larger). The cor- 
responding minimum virial temperature is T^ii^miniz) « 
80, OOO(l-l-z) K for the L22/128 simulation, where the {l + z) 
factor arises from the increasing physical density at higher 
z for fixed virial overdensity. Including the redshift depen- 
dence of the virial overdensity makes the pre-factor slightly 
lower at high redshifts (~ 75000 K) and slightly higher at low 
redshifts (88000 K at 2: = 0) . Any resolved galaxy in the sim- 
ulation resides in a dark matter halo with Tvir > rvir,min(2). 

2.4 Determining Gas Accretion Histories 

At each of our output redshifts (listed below), we use SKID 
to identify all galaxies above our 64mspH baryonic mass 
resolution threshold (6.8 x lO^M© for L22/128). We dis- 
card from consideration any galaxy whose progenitor was 
not resolved at the previous output redshift, since we do not 
want to count simply passing across the resolution thresh- 
old as "accretion." For the remaining galaxies, resolved in 
both outputs, we identify any gas particles that are in the 
galaxy at the later redshift but were not in any resolved 
galaxy at the earlier redshift as smoothly accreted gas. For 
each of these accreted gas particles, we trace back the en- 
tire temperature history since the beginning of the simula- 
tion and record the particle's maximum temperature Tmax. 
These Tmax values will be used to distinguish the cold and 
hot accretion modes. 

To determine the global accretion properties at a given 
redshift, we sum over all the particles that have been 



smoothly accreted between our previous analysis output and 
that redshift. By "smoothly accreted" we simply mean par- 
ticles that were not in resolved galaxies at the previous 
output. For most of our calculations, we include only ac- 
creted gas particles in the statistics. Star particles can also 
be smoothly accreted either because they were in under- 
resolved galaxies at the previous output or because they 
were tidally stripped "field" stars. (Often these "field" stars 
are effectively part of the galactic halo or the intragroup 
or intracluster environment, but SKID still attaches them 
to a galaxy.) However, since all the stars in our simulation 
form within galaxies, the gas that formed these smoothly 
accreted stars was itself accreted onto a (perhaps under- 
resolved) galaxy at a higher redshift, making the Tmax values 
representative of gas accretion at higher redshifts. Including 
accreted star particles in our statistics has only a small im- 
pact, increasing the apparent importance of cold accretion 
at low redshift (see the dotted line in Figure |5| below) . 

Our measured rates of smooth accretion really repre- 
sent the sum of genuinely smooth gas fiow and mergers with 
galaxies below our resolution limit, including systems that 
are not even under-resolved bound objects in the simulation 
but would nevertheless exist in the real universe. In t|4.1l we 
extrapolate the measured distribution of merger mass ratios 
to show that this sub-resolution merging is unlikely to be an 
important correction, so that our measured smooth accre- 
tion rates indeed correspond mainly to smooth gas fiows. 

We use the following redshift spacing for determining 
Tmax values and calculating accretion rates: for < z < 1 
we use Az ~ 0.125, for 1 < z < 3.5 we use ~ 0.25, for 
3.5 < z < 5 we use ~ 0.5, and for higher redshifts we use 
Az ~ 1. In the cosmology we consider here, these redshift 
intervals correspond to time intervals of roughly At = 0.25 
Gyr for z > 3.5, At ~ 0.5 Gyr at z = 1.5 and At ~ 1.5 Gyr 
at 2: = 0. These time intervals are close to the typical infall 
time scales at the corresponding outputs, except at z = 
where they are slightly longer. We have checked the effect 
of using smaller redshift intervals (see H4.2II and find that it 
makes only a small difference to our statistical results. 



3 COLD AND HOT ACCRETION 

Figure ^ shows the distribution of gas particles in the 
temperature-overdensity plane, from the L22/128 simula- 
tion at z = 3. As with previous studies (e.g., Dave ct ^ 
^9^, we identify three important gas phases. The narrow, 
upward sloping locus with p/p < 10 and T < 10^ K consists 
of low density, highly photoionized gas in the intergalactic 
medium (IGM), which is responsible for the Lya forest. The 
tight temperature-density relation in this regime is main- 
tained by the competiti on between adiabatic cooling and 
photoionization heating jHui fc Gnedinlll997l) . The plume 
of particles with p/p ~ 10 - 10* and T ~ lO'^ - lO'^ K is 
comprised of shock heated gas in virialized halos and, at the 
lower density end, around filaments. The narrow, downward 
sloping locus at p/p > 10'^ and T ~ 10* K represents radia- 
tively cooled, dense gas in galaxies. Cooling times at these 
densities are short, so gas remains close to the equilibrium 
temperature where photoionization heating balances radia- 
tive cooling, which is a slowly decreasing function of density. 
According to the conventional sketch described in the 
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Figure 1. Left: Distribution of gas particles in the p — T plane at z = 3, in the L22/128 simulation. One can easily identify three major 
phases: low density, low temperature gas in the photoionized IGM, shock heated overdense gas, and high density, radiatively cooled gas 
within galaxies. Right: Trajectories of 15 particles that accreted onto galaxies shortly before z = 3, illustrating the "cold" (solid lines, 
circles) and "hot" (dashed lines, triangles) accretion modes. Hot mode particles are shock heated above ~ 10^'^ K before cooling, while 
cold mode particles move directly from the diffuse IGM phase to the dense, galactic phase without ever heating above 10^ K. Trajectories 
start at z = 14.9 and end at z = 3. Points mark the individual redshift outputs, which have typical time separations of 0.05-0.1 Gyr. 



introduction, gas that ends up in galaxies starts in the dif- 
fuse phase, enters the shock heated phase, then cools and 
condenses to reach the galactic phase. Figure 03 plots the 
p — T trajectories of 15 randomly selected gas particles ac- 
creted onto galaxies near 2 = 3, starting near z = 15. Some 
of these particles follow just the path described above. How- 
ever, about half of them start in the diffuse IGM phase and 
move directly to the dense galactic phase without ever heat- 
ing above 10^ K. The virial temperature of the smallest re- 
solved halo in the simulation is 3.0 x 10^ K at z = 3, so 
these particles have never been close to the virial tempera- 
ture of any resolved halo. The separation between outputs 
in Figure^ is typically 0.05-0.1 Gyr, and some of the low- 
Tmax particles could shock heat and cool rapidly in between 
two outputs . Howe ver, the cooling radiation arguments of 
iFardal et al.l (1200 ll) show that such "missed cooling" events 
cannot have a large impact. The total energy radiated by 
gas entering the simulated galaxies is of order the acquired 
gravitational potential energy, as expected, and most of this 
energy is radiated by gas with T < 3 x 10'' K. Gas shock 
heated to 10^ K, by contrast, would radiate 90% of its ther- 
mal energy by the time it cools to 10^ K, so it would not have 
much energy left to radiate at low temperature. Counting 
luminosity as Fardal ot al. (2001^) do closes the loophole of 
missing rapidly cooling particles because the high luminosity 
of such particles would compensate for their rarity. 



The two popula tions of trajector ies in Figure^ (also il- 
lustrated in Fig. 4 of lKav et al .^2000t) represent the processes 
we refer to as "cold mode" and "hot mode" gas accretion. 
To quantify the global significance of these two modes, we 
plot in Figure |5| the distribution of Tmax values of accreted 
gas, computed as described in H2.4I Each panel shows the 
Tmax histogram for particles accreted onto resolved galax- 
ies between the previous output and the indicated redshift. 
At high redshifts, the histograms are clearly bimodal: some 
gas enters galaxies after cooling from T ~ 10^ — 10^ K, but 
a large fraction of the accreted gas was never hotter than 
T ~ 10^ K. The bimodal nature of the Tmax histogram and 
the comparable contributions of the two modes to mass ac- 
cretion by galaxies are the primary result of this paper. 

The overall accretion rate declines towards low redshift, 
though a significant amount of growth still occurs because 
more time is available (see Fig. |1| below). The "cold mode" 
region of the histogram broadens towards low redshift, and 
its centroid shifts towards higher Tmax, so by z = the 
distribution is no longer strongly bimodal. As discussed in 
H2.4I we generally exclude accreted stars from our statistics, 
since they reflect physical conditions that prevailed at ear- 
lier times. Dotted histograms in Figure |5| show the effect of 
including accreted stars, which is mainly to extend the low 
redshift histograms to lower Tmax values. 

Figure |3 is similar to Figure |5| but here we plot the 



© 0000 RAS, MNRAS 000, 000-000 



How Do Galaxies Get Their Gas? 7 



CO 

I 

o 







CO 

d 
d 



CO 

d 
d 



CO 

d 

cv 
d 



to 

(JO o 
■Q. 



CO 

d 



d 



.1 I I I I I I I I I I I I I I M I M M I. 




■ Ill h 1 1 1 1 1 1 1 1 1 1 1 1 1 1 



J 1 1 p 1 1 1 [Ti 1 1 1 1 1 1 rji 1 1 L d 




. 1 1 1 1 H 1 1 1 1 1 1 1 1 1 1 1 1 1 rn I . 

z = 0.5- 

.. X 10 ■ 




.1 I I I I I I I I I I I I I I I I I I I I I I. 

z=4.0- 



X 1 



1 1 1 1 1 hT|irrfJTil-i-| M I L 




25-- 



.1 I I I I I I I I I I I I I I I I I I I I I I. 

z = 3.5- 




4567 4567 4567 

log T_ 

Figure 2. Distribution of maximum temperatures of gas accreting onto galaxies. For each particle that was smoothly accreted onto a 
resolved galaxy between the previous output and the plotted redshift, we trace back its history to determine the maximum temperature it 
had at any previous time. Dashed histograms show this distribution in units of /i^MQyr~^Mpc~^ (comoving) per 0.2-dex bin of logTmax- 
Solid histograms have the same shape but are multiplied by an arbitrary constant (as indicated in the panel) to improve visibility. Dotted 
histograms show the effect of including accreted stars in the calculation. 



distribution of Tmax/Tvir, where Tvir is the virial tempera- 
ture of the galaxy's parent dark matter halo, identified as 
described in H2.3I Kt z > 2.5, there is again clear bimodal- 
ity in the distribution, revealing a distinction between gas 
that goes through a strong, virial-type shock before cooling 
and gas that does not. The high temperature portion of the 
histogram peaks at Tmax ~ Tvir at all redshifts. However, 
the spread in Tmax/Tvir values is large, roughly an order- 
of-magnitude, even for this hot mode accretion. The spread 
presumably reflects departures from the spherical geome- 
try adopted in most analytic calculations, plus the effects 
of post-shock adiabatic compression, which can heat gas to 
Tmax > Tvir bcfote cooliug scts iu. The peak of the low tem- 
perature portion of the histogram moves steadily towards 
higher Tmax/rvir with time. In Figure |5| we saw that the 
characteristic Tmax of cold mode accretion also increases 
with time, rising from Tmax ~ 10"-^ K at high redshift to 
Tmax ~ 10^ K at low redshift, but the shift is more rapid 
when we scale to TVir because the virial temperatures of the 
smallest resolved halos are themselves dropping with time 



(see H2.3II . and these halos account for a substantial fraction 
of cold accretion. The clear bimodality of the histograms has 
therefore disappeared by z ~ 1.5 in Figure |3| and at lower 
redshift the cold mode appears as a tail of the distribution 

towards low Tmax/TVir. 

Figures |5| and El show that the cold and hot modes sep- 
arate more cleanly in the distribution of physical temper- 
ature Tmax rather than the scaled temperature Tmax/TVir. 
This result is itself a useful clue to the physics of cold accre- 
tion. If the characteristic temperature of the cold accretion 
mode were determined by gravitationally induced shocks, 
then we would expect it to scale with the characteristic dy- 
namical temperature Tvir. In this case, the minimum in the 
Tmax/T^vir histogram should stay roughly constant, instead 
of rising steadily with time as it does in Figure |3 The min- 
imum of the Tmax histogram, on the other hand, does stay 
nearly constant at Tmax ~ 10^'^ K, near a local minimum 
in the atomic cooling curve where free-free emission takes 
over from helium line cooling (see, e.g.. Fig. 1 of KWH). 
This behaviour could arise if a significant fraction of in- 
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Figure 3. Like Fig.^l but showing the distribution of Tmax/TVir, the ratio of a particle's maximum temperature to the virial temperature 
of the accreting galaxy's dark matter halo. 



falling gas starts at low temperature (T < 10'* K) and gains 
energy through weak shocks or adiabatic compression, in 
which case gas with short cooling times could radiate en- 
ergy without ever heating to high temperature, while gas 
with long cooling times could not. Clearly this is not the 
full story, since the peak of the Tmax histogram does move 
to higher temperature over time, presumably reflecting the 
higher infall velocities associated with larger gravitationally 
induced structures. We will discuss the physics of cold accre- 
tion at greater length in t|6.2l For now, we use the empirical 
evidence of Figure|5|to set the dividing line between the cold 
and hot modes at Tmax ~ 2.5 x 10"^ K, noting that the loca- 
tion of this division is likely a consequence of atomic physics. 
Gravitationally scaled temperatures appear less effective at 
distinguishing hot and cold accretion. 

Figure 121 shows that cold accretion dominates over hot 
accretion at high redshift, then becomes steadily less im- 
portant at lower redshifts. Figure |1] quantifies this depen- 
dence by showing the mean accretion rate in cold and hot 
modes as a function of time and redshift. The left hand 
panel adopts our standard division between cold and hot 
mode, at Tmax ~ 2.5 x 10"" K. The cold accretion rate and 
the total accretion rate rise rapidly together and peak at 



z ~ 3, then decline towards low redshift. Note that this 
Figure shows accretion onto galaxies above our resolution 
threshold; if we went further down the mass function, the 
accretion rate would be flatter at z > 3. The hot accretion 
rate rises later, peaking at 2 ~ 2, and it declines more gently 
thereafter. Cold accretion dominates at z > 3, while hot ac- 
cretion dominates at z < 2. The right hand panel adopts a 
division at Tmax/Tvir = 1. With this deflnition (which is not, 
we think, the most appropriate one), cold accretion domi- 
nates by a large factor at high redshift, and the two modes 
are comparably important at z < 1. 

These global statistics, averaged over the full resolved 
galaxy population, mask an important underlying trend, 
namely a strong dependence of the cold/hot accretion frac- 
tions on galaxy mass. Figure |K| plots the cold accretion frac- 
tion against galaxy baryonic mass at z = 3, 2, 1, and 0. 
Points show the cold fractions for each resolved galaxy, the 
solid curve shows the median value in bins with equal num- 
ber of galaxies (usually ~ 10 — 20 per bin), and the dashed 
curve shows the median hot accretion fraction. The solid and 
dashed curves sum to one by deflnition. There is a strong 
and continuous trend of cold accretion fraction with bary- 
onic mass. The transition mass where cold and hot accretion 
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Figure 4. Redshift history of tlie total smooth gas accretion rate (dashed line), and the rates in cold mode and hot mode (solid and 
dotted lines, respectively). In the left panel, the division between hot and cold modes is at Tmax = 2.5 X 10^ K, while in the right panel 

it is at Tmax/^vir — 1- 



are equally important is nearly constant at Mgai ~ Mq 
for z <2, and slightly higher (Mgai ~ 10^°-^Mq) at 2 = 3. 
We will show in H4.3l that the value of this transition mass is 
insensitive to the numerical resolution of the simulation. The 
overall change of the cold mode accretion fraction with red- 
shift, seen in Figured is largely a consequence of increasing 
galaxy masses at low redshift. 

We obtain similar results if we use, in place of a galaxy's 
baryonic mass, the total mass of its parent halo, as shown 
in Figure |H| At high redshift, the trend with halo mass is 
even cleaner than the trend with galaxy mass; comparing 
Figures 13 and 1^ shows that the low mass galaxies that are 
dominated by hot mode are mostly those that live in high 
mass halos (satellite objects orbiting a larger central galaxy). 
The cold mode fraction in massive halos is slightly larger at 
high redshifts, an effect that is even more pronounced if we 
associate galaxies with FOF halos instead of SO halos. The 
transition halo mass at which cold and hot modes are equally 
important is Afhaio ~ 10^^ * A/0 (slightly higher at z = 3), 
which is the value one would expect given our transition 
baryonic mass of Mgai ~ 1O^°-^M0 if ~ 2/3 of the available 
halo baryons typically end up in the central galaxy. The 
IQ^^ '^Mq transition mass is about a factor of 2 — 3 higher 
than the value found by Birnboim & Dekol (2003) based on 
1-d numerical experiments, a good level of agreement given 
the radically different calculational methods. We will return 
to this point in H6.2I 

Figure |7| characterises the overall contributions of cold 
accretion to the simulated galaxy population by plotting the 
fraction of each galaxy's baryonic mass that was first ac- 
creted in cold mode against the galaxy's present mass. Here 
we consider all particles present in the galaxy, cold gas and 
stars, and trace back each particle's history to determine its 



Tmax value. (Since each star particle comes from a unique 
gas particle, this value is well defined even for stars.) High 
mass galaxies have little or no cold gas accretion (Fig. |3, 
but they are built largely from mergers of lower mass sys- 
tems, and these did accrete a significant fraction of their 
gas in cold mode. Thus, even for the most massive galaxies 
present at 2 = (lower right panel), roughly 40% of the 
baryonic mass was never heated above our adopted thresh- 
old of 2.5 X lO'^K. For galaxies with A/gai < 5 x 1Q^°Mq, 
more than half of the mass was accreted in cold mode, and 
for Afgai < 2 X IO^^Mq, nearly all of it. If we look at higher 
redshifts, the cold mode fraction at fixed mass is slightly 
higher, but the biggest difference is that more of the galax- 
ies are low mass. Every galaxy present in the simulation at 
2 = 3 was built primarily via cold accretion. We conclude 
that, if our simulation predictions are accurate, cold accre- 
tion is a major ingredient in the recipe of galaxy formation. 



4 NUMERICAL ISSUES 
4.1 Sub-resolution merging 

One issue of both physical and numerical interest is whether 
the "smooth" accretion onto our resolved galaxies has a sub- 
stantial contribution from mergers with systems below our 
resolution threshold, since we do not distinguish between 
gas particles that enter individually and those that enter 
in unresolved groups (see H2.41 . If sub-resolution merging 
dominates, then our computed Tmax values could be spu- 
riously low because the gas does its cooling in unresolved 
systems, and the cold accretion process found here might 
simply be analogous to the mergers with low-Tvir objects. 
The best way to address this issue is to run a simulation 
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Figure 5. The cold accretion fraction as a function of galaxy baryonic mass (cold gas + stars), at 2; = 3, 2, 1, and 0. Points show the 
cold fractions of individual galaxies, and solid lines show the median values in bins of baryonic mass. Dashed lines show the median hot 
fraction; solid and dashed curves sum to one by definition. 



that resolves galaxies all the way down to the limit where 
their formation is suppressed by the UV backgrou nd, a 
baryonic mass of about M^ai = W^Mp) ( Efstathiou l |l992: 
Quinn. Katz fc Efstathioul Il996l: IXhoul fc Weinberg IT99a 
Gnedinll200(ir We have one such simulation, L5. 5/128 with 
64 times the mass resolution of L22/128, which has been 
evolved to 2; = 3. We first calculate the smooth accretion 
rates in L5. 5/128 for galaxies above the L22/128 resolution 
limit, counting as smooth accretion all the gas in galaxies 
below this limit just as we do in our analysis of L22/128 
itself. We then calculate the fraction of this inferred smooth 
accretion that is actually due to mergers with galaxies be- 
tween the 1.1 X 10* M0 resolution threshold of L5. 5/128 
and the 6.8 x lO^M© threshold of L22/128. We find that 
this fraction is less than 7% at ^ = 3 and less than 3% 
at z — 4. Our treatment of photoionization in these sim- 
ulations assumes ionization equilibrium (KWH). If we in- 
cluded non-equilibrium heatin g at the epoch of reionization 
jMiralda-Escude fc Ree3ll994ri . then our IGM temperature 
would be somewhat higher, and galaxy form ation would b e 
physically suppressed at a higher mass scale llGnedinll200Cl) . 
further reducing the sub-resolution merger contribution. 

At 2 < 3 we do not have such a high resolution simula- 
tion available, so we must use an indirect method to estimate 
sub-resolution merging. Our approach is similar to that of 
MKHWD, based on extrapolating the mass distribution of 
resolved galaxies that merge with larger resolved galaxies. 
This distribution is plotted in Figure|Hl for the L22/128 sim- 
ulation. We consider all resolved mergers that take place in 



four redshift intervals: z < 0.5, 0.5 <z<l,l<z<2, and 
2 < 2 < 3. We fit a power law dPni/dlogM oc to the 
low mass bins in each panel, where (IPm /d log M is the prob- 
ability per d log M for a galaxy of mass M to merge with a 
larger galaxy. We use the last five bins for the z < 0.5 distri- 
bution and the last four bins for the other distributions. We 
force the fits to go through the measured value in the fourth 
bin, so that we have just one free parameter, a. Best-fit 
slopes range from 0.2 (z < 0.5) to —0.3 {2 < z <3), though 
the small number of bins and substantial Poisson error bars 
make these estimates uncertain. The mass contributed per 
dlogAf scales as M"'^^ , so high mass satellites dominate 
merger growth if a > —1, as we find at all redshifts. 

Integrating our power law fits from the resolution limit 
Mgai = 6.8 X 10^ Mq down to Afgai = 10* Mq yields an up- 
per limit to the amount of sub-resolution merging, since the 
merger mass distributions appear to be turning over steadily 
and are likely to fall below our extrapolations. With this es- 
timate, the ratio of mass added to the galaxy population by 
sub-resolution mergers to mass added by resolved mergers 
is approximately 0.1 — 0.15 at z = — 0.5, 0.20 — 0.35 at 
z = 0.5-1, 0.25-0.45 at z = 1-2, and 0.5-2.0 at z = 2-3. 
The quoted ranges reflect the la statistical uncertainties in 
the fitted slopes. Sub-resolution merging could thus make a 
significant contribution to merger rates at higher redshifts, 
but because merger rates are much lower than smooth ac- 
cretion rates (see Figure lT^ belowl. it still does not substan- 
tially change the inferred accretion rates. The exception is 
at 2; = 2 — 3, where the maximal correction above could raise 
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Figure 8. The mass spectrum of galaxies that merge with larger 
galaxies, in four redshift intervals as indicated. We fit power laws 
to the low ends of these distributions as described in the text. 
Straight lines show these fits, which have logarithmic slopes a as 
listed. 



the merger rate from 25% of the smooth accretion rate to 
150% of the smooth accretion rate. In other words, up to half 
of our estimated smooth accretion rate at 2: = 2 — 3 could be 
due to sub-resolution merging. However, at 2: = 3 we have 
direct evidence from the comparison of the L5.5/128 and 
L22/128 simulations that the unresolved merger rate is only 
~ 7% of the smooth accretion rate, implying that our power 
law extrapolation greatly overestimates the sub-resolution 
merging at high redshift, where the typical galaxy mass is 
only a few times the resolution limit. We hope to eventually 
redo our analysis with higher resolution simulations that 
model larger volumes than L5. 5/128 and continue all the 
way to z = 0, but for now we conclude that the quantity 
we estimate as the smooth accretion rate in the L22/128 
simulation does represent truly smooth gas accretion, with 
only a modest contribution from sub-resolution mergers. We 
thus confirm the conclusion of MKHWD about the mecha- 
nisms of galaxy growth: typical galaxies gain most of their 
mass through smooth accretion of gas rather than mergers 
with smaller galaxies.^ Furthermore, the cold accretion pro- 

This agreement is, however, somewhat fortuitous. MKHWD ap- 
plied the power law extrapolation method to the L50/144nb simu- 
lation, but their Figure 9 incorrectly labels the merger mass distri- 
bution as being per dM instead of d log M (the quantity they actu- 
ally calculated). This translates to an error in the power law coeffi- 
cient, which caused MKHWD to underestimate the sub-resolution 
merging. If MKHWD had correctly applied their power law fit to 
their data, they would have concluded that sub-resolution merg- 
ing exceeds the true smooth accretion rates. The higher resolution 
of the L22/128 simulation allows us to sec a turnover in the merger 
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Figure 7. Contribution of cold and hot accretion to total galaxy masses at z = 3, 2, 1, and 0. Points show the cold mode fraction of 
each galaxy, i.e. the fraction of its mass with Tmax < 2.5 X 10^ K, as a function of baryonic mass. The solid line is the median cold mode 
fraction, while the dashed line is the median hot mode fraction; by definition, the two fractions sum to one. 



cess that is the focus of this paper is not a consequence of 
galaxies accreting gas that has cooled in halos of low virial 
temperature. 

4.2 Time resolution 

A second technical issue is the time interval between our 
analysis outputs. Since Tmax is by definition the maximum 
temperature of a particle during its pre-accretion history, 
our analysis at discrete outputs necessarily yields a lower 
li mit to the tru e Tmax value of any particle. The arguments 
of iFardal et alJ (1200 ll) imply that we cannot be drastically 
underestimating Tmax by missing rapidly cooling particles, 
since these would then have high luminosities that would 
cause the cooling radiation to emerge mainly at high tem- 
peratures, which it does not. However, we have also carried 
out a direct check by repeating our analysis at several z > 1 
outputs using our standard time interval to determine which 
particles are smoothly accreted but determining Tmax val- 
ues at smaller intervals (by a factor ~4atz = 3,~8at 
2 ~ 2, and ~ 20 at z = 1). With this reanalysis, we find 
that 15-25% of the smoothly accreted gas that was previ- 
ously attributed to cold mode accretion is actually hot mode 

mass distribution, which is why we find a small sub-resolution 
contribution from correct application of the same method. 



accretion, with the larger fraction applying at low redshift 
when the cold fraction is itself small. Further investigation 
shows that most of the "missed" hot particles shock heat 
well inside the virial radius, where the density of the gas 
is high and the cooling times are short, so they are physi- 
cally distinct from the "classic" hot mode particles heated at 
the virial shock. With higher time resolution, the Tmax dis- 
tribution of hot mode particles also shifts towards slightly 
higher temperatures. Note, however, that if we use higher 
time resolution for all aspects of our analysis, rather than 
just for computing Tmax values, then we can actually get an 
increase in the cold mode fraction at higher redshift, since 
the calculation now incorporates low mass galaxies that were 
previously omitted because they were below the resolution 
limit at the preceding redshift output. Using the finer time 
resolution described above to identify galaxies and compute 
Tmax values yields a net increase in the cold mode fraction 
of 30% at z ~ 1, no significant change at z ~ 3, and a small 
drop at higher redshifts. 

We conclude that the finite time resolution in our anal- 
ysis may cause errors of ~ 0.1 — 0.2 in our estimates of global 
cold and hot accretion fractions. This weak dependence on 
time resolution makes sense if we ask how far inside a halo a 
particle needs to be to get shock heated and cool before the 
next output redshift, and thus be misidentified as cold ac- 
cretion. Our regular analysis interval is At ~ 1.0 Gyr at low 
redshifts and At ~ 0.2 Gyr at high redshifts. Table 2 lists 
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348 


69.6 
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294 
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150 


75.2 


15 
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Table 2. Cooling times in Myr of gas at specified overdensity and redshift, for starting temperatures T = 10^ K (top) and T = 2 X 10^ K 
(bottom). 



the time required for gas to cool from 10^ K or 2 x 10® K 
(typical halo virial temperatures) down to 10^ K, at various 
redshifts and overdensities. We compute these times by in- 
tegrating over the actual cooling function. For our standard 
time interval, a particle heated to 2 x 10® K cannot cool to 
T = 10^ K between two redshift outputs unless its overden- 
sity is greater than 500, corresponding to less than 1/3 of the 
virial radius for an isothermal halo. For a particle heated to 
T = 10^' K the cooling time at overdensity 500 is comparable 
to our standard time interval, which means that the cooling 
time at the virial radius is much longer than the time interval 
between our simulation outputs. For our finer time interval, 
a shock heated particle would only cool between two out- 
puts for an overdensity of several thousand, requiring it to 
be very close to the galaxy's cold baryon component. 

4.3 Numerical Convergence 

The most serious numerical concern is that our basic result, 
the existence of a cold accretion mode that makes a major 
contribution to galaxy formation, is somehow an artifact of 
the simulations' finite mass resolution. For example, numer- 
ical broadening of shocks might allow gas to radiate energy 
while it is passing thr ough the shock, inst ead of first heat- 
ing and then cooling jHutchings fc Thomas. .2000,') . We can 
address this issue by comparing simulations with different 
mass resolution to check the stability of the results. 

The top panels of Figure |U] compare the Tmax distribu- 
tions of gas accreted at z = 3 and z — 4 in the L22/128 
simulation and the higher resolution, Lll/128 simulation. 
In Lll/128, we only consider the accretion onto galaxies 
above the 6.8 x 10^ M© resolution limit of L22/128. These 
objects are resolved by at least 512 gas and star particles in 
Lll/128. Despite the factor of eight difference in mass res- 
olution, the two histograms agree nearly perfectly at z = 4, 
and they show only modest differences at 2: = 3. Since the 
simulation volumes are also different, we expect some dif- 
ferences between these simulations just because of the dif- 
ferent structures they contain. The lower panels of Figure 
compare the Lll/128 and Lll/64 simulations, which have 
identical initial conditions (except for the additional high 
frequency modes in Lll/128) and difi'er only in mass resolu- 
tion. The change of simulation volume slightly improves the 
agreement at z = 3, but it appears that most of the differ- 



ence between L22/128 and Lll/128 is an efi'ect of resolution 
rather than simulation volume. Higher resolution lowers the 
overall gas accretion rate by ~ 30%. Some of the gas ac- 
creted in the lower resolution simulation is converted into 
stars within galaxies below the resolution limit of the lower 
resolution simulation, but that are resolved in the higher 
resolution simulation, before it is accreted onto a galaxy 
above the 512-particle threshold. This can explain a signifi- 
cant part of the differences in the cold mode accretion rates. 
However, numerical resolution effects on the gas cooling are 
probably responsible for most of the differences in the hot 
mode accretion rates and the remaining cold mode differ- 
ences. These numerical effects must depend on the galactic 
mass, since at z = 4, where galaxies are significantly smaller 
on average, the differences are smaller. Since the difference 
between L22/128 and Lll/128 is fairly small and mainly at- 
tributable to mass resolution, we conclude that the volume 
of the L22/128 simulation is large enough to give reasonably 
accurate statistics on the cold and hot accretion fractions, 
at least at z = 3. 

Our most powerful test of numerical convergence comes 
from examining the galaxy mass dependence of the cold ac- 
cretion fraction, as shown for the L22/128 simulation in Fig- 
ure 13 Figure llOb repeats the 2 = 3 panel of Figure |3 and 
Figure ITOb shows the same result for the Lll/128 simula- 
tion. While the higher resolution simulation contains few 
high mass galaxies because of its small volume, the tran- 
sition between cold and hot mode domination occurs at 
A/gai = 1O^°-^M0, close to the transition in L22/128. Our 
highest resolution simulation, L5. 5/128, contains only nine 
galaxies above the L22/128 mass resolution limit at 2 = 3 
(and only one at 2 = 4), so a comparison of similar form 
is not very useful. However, we can check the mass depen- 
dence of the cold accretion fraction by following the largest 
(and thus best resolved) galaxy in the simulation volume. 
Figure llUfc plots the cold and hot accretion fractions of 
this largest galaxy as it grows from Mgai = 1.4 x 10^'' M0 
(8000mspH) at 2 = 6 to Mg^i = 7.4x lO^^M© (43, OOOmspu) 
at 2 = 3, the redshift of our last simulation output. The tran- 
sition between hot and cold mode dominance occurs close to 
the 2 = 4 output, when the galaxy mass is Mgai = 10^'^'^ Mq. 
This agrees well with the transition mass at 2 = 3 found in 
Lll/128 and L22/128 (and agrees even better with the 2 = 4 
transition mass, which is 10^°-'^^ Mq for L22/128). The sec- 
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Figure 9. Dependence of the Tmax distribution on numerical resolution. Upper panels compare results from our standard, L22/128 
simulation (solid histogram) to results from the higher resolution, Lll/128 simulation (dotted histogram) at 2 = 3 (left) and 2 = 4 
(right). Differences can be partly attributed to resolution and partly to simulation volume. Lower panels isolate the resolution effect by 
comparing simulations Lll/128 and Lll/64, which have the same initial conditions and volume and differ only in mass resolution. 



ond most massive galaxy in L5.5/128 does not reach the 
cold/hot transition mass until 2 = 3, the final output, so it 
is not useful for this test. The agreement in transition mass 
over simulations that span a factor of 64 in mass resolution is 
the strongest single piece of evidence that the importance of 
cold accretion in our simulations is not an artifact of limited 
resolution. 



4.4 Influence of the UV Background 

Figure lllb illustrates the impact of the UV background 
on Tmax distributions, comparing results from Lll/64 and 
Lll/64nb at 2 = 3. The no-background simulation shows 
a somewhat reduced hot mode, and the cold accretion is 
sharply peaked at Tmax ~ 10* K. The UV background is the 
only difference between these two simulations. The different 
amplitude and shape of the histogram at low Tmax is there- 
fore driven by the difference in input physics, most likely 
the fact that photoionization reduces the strength of the 
coUisional line peak in the cooling curve at moderate over- 
density (see KWH, Fig. 2). The transition mass at which 



cold and hot accretion rates are equal is consistently higher 
in simulations with no UV background, increasing the global 
contribution of cold accretion. 

L22/128 is the highest resolution simulation we have 
evolved to z = 0. All of our lower resolution simulations 
have been evolved without a UV background, because pho- 
toionization at lo w mass resolution spuriously supp resses 
galaxy formation ijWeinberg. Hernauist. fc Katjll997f l. Fig- 
ure lllb compares the z = Q Tmax histograms from L22/128, 
L22/64nb, and L50/144nb. We now consider only the ac- 
cretion onto galaxies above the 64-particle threshold of the 
lower resolution simulations, Mgai = 5.4 x 10^" M©, which is 
higher than the cold/hot transition mass of ~ 3 x 10^" 
found for L22/128 at z = (see Fig. O. As a result of this 
higher threshold mass, the low temperature end of the Tmax 
histogram for L22/128 is strongly suppressed relative to the 
corresponding histogram in Figure |2| The no-background 
simulations, by contrast, show a significant amount of cold 
accretion, with Tmax values peaked just above 10* K. Al- 
though the simulations differ in both numerical parameters 
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Figure 10. Resolution tests of the cold accretion fraction, (a) 
Cold fraction vs. galaxy baryon mass for the L22/128 simulation 
at z = 3, repeated from Fig. |S] Points show individual galaxies, 
the solid line shows median cold fraction, and the dashed line 
shows median hot fraction, (b) Same plot for the Lll/128 simu- 
lation, at z = 3. (c) Cold fraction vs. galaxy baryonic mass for 
the single largest galaxy in the L5. 5/128 simulation at redshifts 
2 = 5.5, 5, 4.5, 4, 3.5, and 3 (each point is a redshift, with the 
highest redshift corresponding to the lowest mass). Dashed line 
shows the hot fraction. While the mass resolution of this sim- 
ulation is 64 times higher than that of L22/128, the crossover 
between cold mode domination and hot mode domination occurs 
at nearly the same baryonic mass. 



and input physics, the similarity to Figure [TTh implies that 
the UV background is the cause of this difference. 

The L22/64nb and L50/144nb simulations also show 
significantly more hot accretion at z = than L22/128; 
a similar effect is seen at z = 3 (not plotted). Here we sus- 
pect that the coarse resolution of the no-background simu- 
lations is to blame, since photoionization has minimal effect 
on cooling at high temperatures and the differences are sim- 
ilar to those seen in Figures andEJ; (at 2 = 3). Finally, we 
note that the reasonably good agreement between L22/64nb 
and L50/144nb suggests that the 22/i~^ Mpc box is large 



enough to give representative statistics for the cold/hot frac- 
tion even at z = 0, except perhaps at very large Tmax. 

4.5 Other simulations 

Within our own set of simulations, we find that the exis- 
tence of an important cold accretion mode is insensitive to 
mass resolution or simulation volume, though the specifics of 
the Tmax distribution are significantly affected by the pres- 
ence of a UV background. SPH simulations of the formation 
of individual galaxies achieve still higher resolution, with 
thousands of particles in the central o bject, and they also 
show a large fraction of cold accretion IIKatz &: Gunnlll99ll : 
Abadi et al. 2003r). Further evidence of the robustness of the 
result comes from simulations of volumes like those consid- 
ered here p erfor med with entirely independent SPH codes. 
iKav et all i200Ql . studying a 2 x 32^ particle simulation 
of a 10h~^ Mpc box performed by the Virgo consortium, 
find that only ~ 11% of the gas accreted by their simu- 
lated galaxies by z = was ever hotter than 10"" K. We 
have recently analysed a simulation, kindly provided by L. 
Hernquist and V. Springel, that uses Springel & Hernquist's 
(2002) entropy conserving formulation of SPH, which makes 
the numerical treatment of shock heating significantly dif- 
ferent from that in our code . This simulation, run Q3 of 
ISprineel fc Hernauisd ll2003bft . represents a Wh ^ Mpc box 
with 2 X 144'^ particles, and is evolved to z = 2.75. The sim- 
ulation also includes strong f eedback from stellar w i nds, us - 
ing the methods described bv lSpringel fc Hernauist| (l2003a^ . 
Despite the se differences in numerical m ethod and input 
physics, the ISpringel fc HernauistI ll2003tJ) simulation pro- 
duces a Tmax histogram at high redshifts that closely resem- 
bles those shown in Figure |5| with a large fraction of gas 
accretion occurring at low temperatures. The transition be- 
tween cold mode dominance and hot mode dominance occurs 
at a halo mass similar to that found in Figure |S| However, 
the strong feedback makes galax y baryonic masse s signifi- 
cantly smaller in the iSpringel fc Hernguisti ^l2003b^ simula- 
tion than in our simulation (numerical differences may also 
play a role), so the baryonic mass at which this transition 
occurs is lower than that in Figure |S] by a factor ~ 4. In 
sum, galaxy masses and the amount of accretion are lower 
in the lSpringel fc Hernauisli (i2003b) simulation, but the ex- 
istence of a cold accretion mode and the relative amounts of 
cold and hot accretion at a given halo mass appear entirely 
consistent with our results. 

These internal and external comparisons strongly sug- 
gest that the bimodal Tnax histogram is a genuine physi- 
cal result, at least given the physical assumptions of these 
simulations. Since the results quoted above are all based 
on SPH simulations, albeit with different numerical imple- 
mentations, it would be desirable to confirm the result with 
Eulerian grid or adaptive mesh hydro simulations, which 
typically use shock capturing metho ds instead of ar tificial 
viscosity to implement shock heating. ' Cen fc Ostrikeil a999. 
Fig. 4) find a broad temperature distribution for cooling ra- 
diation in their Eulerian hydro simulation, suggesting that 
much of the gas in this simulation also cools without ever 
reaching T > 10^ K. A. Kravstov (private communication), 
in adaptive mesh simulations of individual forming galaxies, 
finds that much of the gas penetrates along cold filaments far 
inside the virial radius without being shock heated. These 
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Figure 11. Influence of numerical resolution, simulation volume, and the UV background on the Tmax distribution, (a) Comparison 
of Lll/64 and Lll/64nb, which differ only in the presence of a UV background, at z = 3; (d) Volume effect comparisons of L22/128, 
L22/64nb and L50/144nb at 2 = 0. 
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Figure 12. Star formation rate per unit comoving volume (solid 
line) for all resolved galaxies in the simulation compared to 
smooth gas accretion (dashed line) and merger accretion rates 
(dotted line). The merger mass accretion rates include total mass 
gain in merger events (gas and stars). 



simulations cannot presently measure the Tmax histogram 
itself because they do not track gas from cell to cell and 
therefore do not record the history of gas that ends up in 
galaxies, but this measurement can be carried out in the fu- 
ture by adding a population of "tracer" particles that foUows 
the gas flow (Kravtsov, private communication). 



5 THE ENVIRONMENT DEPENDENCE OF 
ACCRETION AND STAR FORMATION 
RATE 

Figure 1121 plots the globally averaged star formation rate 
(SFR) and the globally averaged contributions of smooth 
accretion and mergers to galactic mass growth, as functions 
of time in the L22/128 simulation. At high redshifts the 
smooth accretion rate dominates the merging growth rate 
by a large factor, while at z < 1 these two processes have 
comparable global rates. Of course, galaxies that gain most 
of their mass via mergers do exist in the simulation, but 
the growth of the typical galaxy is dominated by smooth 
accretion. MKHWD found that the global star formation 
history in the L50/144nb simulation closely tracks the gas 
accretion history and does not track the galaxy merger rate. 
Figure 1121 confirms this conclusion using our higher reso- 
lution, L22/128 simulation — the correlation between the 
accretion rate and SFR curves is even tighter than the one 
found by MKHWD. The two curves are offset by ~ 0.5 Gyr, 
reflecting the average time that it takes gas to be converted 
into stars after it is accreted by a galaxy. Ho wever, with 
our observationally motivated assumption (e.g., iKennictitl 
1998) that the star formation rate is an increasing function 
of gas density, any steadily accreting galaxy will form stars, 
and mergers are not needed as triggering events. Indeed, the 
global star formation rate must roughly track the global ac- 
cretion rate for any star formation model in which galaxies 
do not build up large, quiescent gas supplies. Mergers could 
play a role in causing rapid bursts of star formation by driv- 
ing galactic gas to higher densities, and our simulations may 
underestimate this effect because of their limited resolution. 
However, such merger- induced acceleration would have lit- 
tle effect on the globally averaged SFR because the merging 
galaxies simply consume the same amount of gas on a some- 
what shorter timescale. 

If our input physics is approximately correct, then 
understanding the redshift and environmental dependence 
of the cold and hot accretion rates is tantamount to un- 
derstanding the redshift and environmental dependence of 
galactic scale star formation. In this section we focus on the 
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environmental dependence; we return to the redshift depen- 
dence in § 16.31 and § 16.41 Figure 1131 plots the mean SFR 
and accretion rate per galaxy as a function of environment, 
specifically the comoving space density of resolved galaxies 
computed using SPH-style kernel averaging on a galaxy-by- 
galaxy basis with a variable size spli ne kernel whose size i s 
chosen to enclose 15 neighbours (as in lWeinbere et al.l2no4l . 
The SFR and accretion rate track each other as a function 
of density, just as the global rates track each other as a func- 
tion of redshift. At z = 3, the average accretion rate is nearly 
independent of environment, but at lower redshifts there is 
a break towards lower accretion rates and star formation 
rates at ngai ~ 1 /i"^ Mpc""^. Points in Figure IT!T1 show the 
median accretion rate, computed in bins containing equal 
numbers of galaxies (while the mean is computed in bins of 
constant A log n = 0.2). Comparing the medians and means, 
one can see that the accretion rate in high density regions 
is dominated by the few galaxies with large accretion rates; 
most galaxies in these regions have very low accretion and 
correspondingly low star formation. 

The dotted and dashed lines in Figure HTH show the sep- 
arate contributions of cold and hot accretion to the mean ac- 
cretion rate. There is a clear trend, increasingly pronounced 
at lower redshift, of cold mode domination in low density 
environments and hot mode domination in high density en- 
vironments. The transition between the two modes shifts 
to progressively lower comoving densities at lower redshifts. 
Cold mode dominates in all but the highest density en- 
vironments at z ~ 3, but at z = it dominates only 
in the lowest density environments and is insignificant at 
Wgai ^ 1 h'^ Mpc""^. These trends can be largely explained 
by the galaxy mass dependence of the cold/hot accretion 
fraction: low density regions are populated by lower mass 
galaxies, which tend to have more cold accretion. However, 
when we examine the trends separately for low and high 
mass galaxies, we find that there is also a direct environmen- 
tal effect, in that low mass galaxies are often hot mode dom- 
inated in high density environments, while they are always 
cold mode dominated in low density environments. This dif- 
ference reflects the higher characteristic halo masses in dense 
environments; as noted in our discussion of Figures |S| and |S1 
low mass galaxies in high mass halos tend to be hot mode 
dominated. The rare high mass galaxies in low density en- 
vironments are still hot mode dominated, but they tend to 
have somewhat higher cold accretion fractions than those in 
denser regions. 

We can c ompar e our z = predictions to the results of 
Gomez et al. 1120031 hereafter G03), who use Ha equivalent 
widths to infer the star formation rates of a volume limited 
sample of galaxies from the SDSS. To mimic their selection 
of galaxies above an absolute magnitude limit M* + 1, we 
choose galaxies above a baryonic mass threshold that yields 
the same galaxy number density, roughly the 120 most mas- 
sive galaxies in our 22.222/1"^ Mpc box. Like G03, we cal- 
culate the surface density around each galaxy by finding 
the distance rio to the tenth nearest neighbour in projec- 
tion, S = 10/(7rrio), after eliminating galaxies further than 
±1000 km s~^ in redshift (making use of our periodic bound- 
ary). To improve our statistics, we combine the results from 
the three orthogonal projections of the box. 

Figure lT^ plots the median star formation rate and the 
25th and 75th percentiles of the SFR distribution in bins of 



galaxy surface density containing equal number of galaxies, 
similar to Figure 4 of G03. The simulation reproduces the 
principal qualitative result of G03, namely a break in the 
SFR vs. E relation at E ~ 1 h^^ Mpc~^, most prominent 
for the 75th-percentile curve, with a slowly increasing SFR 
below the break and a rapidly declining SFR at higher sur- 
face densities. G03 also examine the dependence of SFR on 
the distance to the centre of the nearest group or cluster. 
They emphasise that the break in SFR occurs at 3 — 4 times 
the cluster virial radius, implying that it is not an effect as- 
sociated with virialization or the intracluster medium. We 
have just one prominent cluster in the L22/128 simulation, 
with virial mass 2.7x lO^*Af0. Figuredt plots SFR against 
projected distance to the centre of this cluster (defined by 
the position of the most bound particle), in units of the clus- 
ter virial radius _Rvir = l-lh~^ Mpc, for comparison to GOS's 
Figure 6. The breaks in the simulated SFR curves also occur 
beyond the cluster virial radius, at roughly 2_Rvir. At least 
part of the drop at large radii is caused by the presence of a 
large galactic group at distances between 1 and 2 -Rvir, which 
dominates the statistics. A more detailed comparison would 
require a larger simulation volume that would allow us to 
mimic G03's group selection method more carefully. Since 
the SFR tracks the gas accretion rate in our simulations, a 
drop in SFR within high density regions at low redshifts is a 
direct consequence of a diminished gas supply and does not 
require cluster specific processes that either rapidly consume 
the existing gas (e.g. galaxy harassment) or remove it from 
galaxies (e.g., ram pressure stripping). 

G03 note that their SFR values should be multiplied 
by a factor ~ 5 to account for the 3" aperture of the SDSS 
fibers, but even so they are significantly lower than the rates 
predicted by the simulation. For example, the plateau of the 
75th-percentile SFR curve is at ~ 2 Mq yr"^ gar^ in G03 
(including the factor of five) and ~ 5 AIq yr~^ gal^^ in 
the simulation. The 25th-percentile curve is consistent with 
zero SFR at all E in G03, while the simulation predicts 25th- 
percentile rates ~ 1.5 Mq yr~^ gal~^ in low density environ- 
ments. G03 include negative estimates of the SFR derived 
from Ha flux in their median calculation, which slightly low- 
ers their SFR values, but probably not by a large factor. 
Limited mass resolution of our simulation may have some 
effect in making the predicted SFR and accretion rates ar- 
tificially high (see § 14.311 . In any event, we are encouraged 
that the simulation reproduces one of the most striking fea- 
tures of the observed correlation between SFR and envi- 
ronment, a break at surface densities ~ 1 /175 Mpc~'^ that 
lie beyond the virial radii of rich galaxy gro ups, a result 
found independently by G03 in the SDSS and bv lLewis et alJ 
J2002D in the 2dF Galaxy Redshift Survey. (The SDSS re- 
sults oflKauffmann ct al. 2004] also appear consistent with 
this finding, although they do not present them in this form.) 
This break occurs in a density regime where the accretion 
rates are dominated by the hot mode (see Fig. I13II . and 
galaxies above G03's M* + 1 threshold are mostly hot mode 
dominated in any case. Thus, although these observations 
test some aspects of the simulation predictions, they do not 
test the existence of an important cold accretion mode or the 
transition of the dominant accretion mode from cold to hot 
with increasing galaxy density. Such effects should become 
more apparent if the galaxy sample is extended to include 
lower luminosity systems. 
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Figure 13. Dependence of accretion rates and star formation rates on local galaxy number density (comoving), at 2; = 3, 2, 1, and 
0. Solid lines show the mean total accretion rates in bins of Alogn = 0.2, and dotted and dashed lines show the mean cold and hot 
contributions, respectively. Points show the median total accretion rate in bins containing equal galaxy numbers. Dot-dashed lines show 
mean star formation rates, which track the total accretion rate. 




Figure 14. (a) Star formation rate vs. galaxy surface density at z = 0. Curves show the 25th-percentile, median, and 75th-percentile of 
the SFR per galaxy in bins of surface density, (b) Same as (a), but galaxies are binned as a function of distance from the centre of the 
largest cluster in the L22/128 simulation, in units of the cluster's virial radius. 
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6 DISCUSSION 

6.1 Physics of the Hot Mode 

In our simulations, liot accretion (Tmax > 2.5 x 10^ K) 
dominates over cold accretion in galaxies of baryonic mass 
Mgai > W^^-^Mq (Fig.El or in halos with Mhaio > lO^^M© 
(Fig. ISJ. Roughly half of the total baryonic mass of high 
mass galaxies (Mgai > IO^^Mq) originates from hot accre- 
tion, while lower mass galaxies are built predominantly 
by cold accretion (Fig. |7|l. The hot accretion mode in 
the simulations corresponds reasonably well to the accre- 
tion envisioned in the standard picture of galaxy forma- 
tion that underlies most semi-analytic models: gas shock 
heats to roughly the halo virial t emperature before cooling 
and settling into a galaxv iRees fc Ostriker 1977 : Silk 19771: 
IWhite fc Ree3ll97d lwhite fc Frenklll99ll) . 

The actual distribution of Tmax/Tvir for hot mode ac- 
cretion is broad, roughly an order of magnitude (Fig. |3Ji. 
This broad distribution reflects a combination of depar- 
tures from spherical symmetry, hierarchical assembly, and 
post-virialization heating. With aspherical infall, gas can 
hit an accretion shock with a range of velocities and thus 
heat to different temperatures. Some gas may be shocked to 
Tmax < 2vir in filamentary structures beyond the virial ra- 
dius, then enter the halo without experiencing a true virial 
shock. When an individual galajcy falls into a group or clus- 
ter, some of its halo gas is heated to the higher virial tem- 
perature of the new halo and escapes the galaxy potential 
well, but other gas remains at the galaxy's virial temper- 
ature and may be accreted before it is heated or stripped. 
Since we calculate the ratio of Tmax to the virial temper- 
ature of a galaxy's host halo at the time of accretion, hot 
mode gas accreted in this way will also have low Tmax/Tvir. 
Finally, gas that heats by an accretion shock near the virial 
radius can be heated further by adiabatic compression as 
it flows inwards and increases its density. We find that a 
significant fraction of the hot mode gas has T significantly 
below Tvir (but still much higher than the typical cold mode 
temperature) in the outer regions of the halo and does not 
reach Tvir until it is inside ~ 0.5-Rvir. 

Figure 1151 shows the same p — T trajectories as Fig- 
ure0D, but now the horizontal axis represents physical den- 
sity (scaled to the mean baryon density at z = 3) rather 
than overdensity. In this plane, adiabatic compression moves 
a particle along a line logT oc (2/3) log p. In the diffuse 
IGM phase, particles move to lower densities and tempera- 
tures as the universe expands, along a locus of slope close to 
0.6. Hot mode particles then experience sudden shock heat- 
ing, which is often followed by an extended phase of heating 
along a nearly adiabatic trajectory, before eventually cooling 
to 10* K and high overdensity. The fact that shock heating 
is spread over several outputs (which are typically spaced 
by ~ 0.05 — 0.1 Gyr) is a consequence of the flnite numerical 
width of the shocks. We will discuss the behaviour of cold 
mode particles in this plot below. 

We have also examined the distribution of cooling times 
in halos at various redshifts. If we calculate the amount of 
hot gas that would cool between two outputs in the absence 
of further heat input, we get a value that exceeds the actual 
hot mode accretion rate by a factor ~ 2. Therefore, adiabatic 
compression or further shock heating must be sufficient to 
prevent roughly half of this gas from cooling. 
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Figure 15. Density-temperature trajectories for the same parti- 
cles shown in Fig. ^3, but with the horizontal axis representing 
physical density (scaled to the cosmic mean density at 2 = 3) 
instead of overdensity. The line segment at the bottom has the 
slope dlogT/dlog p = 2/3 expected for pure adiabatic evolution. 

6.2 Physics of the Cold Mode 

Why does a substantial fraction of gas enter galaxies without 
ever heating close to the virial temperature of the galaxy's 
parent dark matter halo? Our results do not provide a com- 
prehensive answer to this question, but they do offer a num- 
ber of clues that allow us to build a physical picture of how 
cold accretion occurs. As discussed in the division be- 
tween hot and cold modes occurs at roughly constant Tmax 
rather than constant Tmax/Tvir, which suggests that it is 
determined mainly by cooling physics rather than by shock 
heating. The strong dependence of the Tmax distribution for 
cold mode gas on the presence of a UV background (Fig. llH 
further supports this inference. The ratio of cold to hot ac- 
cretion is tightly correlated with galaxy mass or halo mass, 
and the transition between cold mode dominance and hot 
mode dominance occurs close to (though somewhat above) 
th e halo mass predicted by t he spherically symmetric models 
of lBirnboim fc Dekell (l2003h . The ratio of cold to hot accre- 
tion also depends on the density of the environment, at least 
at low redshift, and although this dependence largely reflects 
the relative galaxy mass functions in high and low density 
regions, low mass galaxies are dominated by hot accretion 
in high density environments and by cold accretion in low 
density environments. 

The trajectories of particles in the p — T plane offer 
further insight into the physics of cold accretion. Figure^ 
shows that typical cold mode particles, after an initial phase 
of adiabatic cooling, experience most of their heating at 
overdensities ~ 3 — 30, then steadily increase their overden- 
sity while cooling from slightly above 10* K to 10* K. Fig- 
ure^] plotted in physical density rather than overdensity, 
shows that the heating phase of these trajectories is very 
steep, implying that it is produced by shocks not by adia- 
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batic compression. The characteristic shock velocities must 
be w ~ {kT / ^.m^f^ ~ 10 - SSkms'^ (T ~ 10* - 10^ K) 
and are probably around filaments. 

To better understand the origin of the hot /cold distinc- 
tion, it is useful to look for differences in the physical state 
of gas before it is accreted in one of these two modes. Fig- 
ure llfcil plots the average overdensity and temperature his- 
tories of gas that is accreted onto galaxies a.t z = 2. Solid 
lines show histories for cold mode particles and dashed lines 
for hot mode particles. Averaging over many particles tends 
to smooth the density and temperature evolution, since in- 
dividual particles experience their rapid heating phases at 
different times. The important result in Figure llfcil is that 
cold mode particles initially have higher overdensities than 
hot mode particles, at least on average. Because of the cor- 
relation between density and temperature in the photoion- 
ized IGM (Fig.0, cold mode particles also start at slightly 
higher temperature. The differences are fairly small, and be- 
come somewhat larger if we focus only on the lower-Tmax 
cold mode particles. The general appearance of Figure [T^ is 
similar at other redshifts. 

Since overdensities ~ 3 — 30 are associated with filamen- 
tary structures, while halo overdensities are > 100, these 
p — T trajectories suggest that cold mode particles experi- 
ence much of their heating in filamentary structures, while 
hot mode particles are heated at halo virial shocks. Fig- 
ure ll7l provides anecdotal evidence to support this idea. The 
left panels show an example of a cold mode dominated halo 
(Mhaio = 2.6 X 10^^ Mq) containing a well resolved galaxy 
(Mgai = 1.6 X 10^° Mq « 10,000mspH) in the L5.5/128 
simulation at z — 5.52. The right panels show this same 
halo and its surroundings at z — 3.24, when it is well into 
the hot mode dominated regime (Mhaio = 1.26 X 10'^ Mq, 
Afgai = 7.35 X 10^° Mq). Particles are colour coded by over- 
density in the upper panels and by temperature in the mid- 
dle panels, all of which are 4i?vir on a side and 2_Rvir thick. 
The lower panels show a zoomed region of the temperature 
plot, 1-Rvir on a side and 0.5i?vir thick. Lines attached to 
particles show their projected velocities. In the left hand 
panels, green particles represent cold mode gas that will be 
accreted by the central galaxy by z = 5, 130 Myr later. In 
the right hand panels, they show hot mode gas that will be 
accreted hy z = 3, 190 Myr later. 

The z — 3.24 picture corresponds reasonably well to the 
conventional notion of galactic gas accretion, with accret- 
ing material coming from the inner regions of a hot, quasi- 
spherical halo. However, even at this redshift cold filaments 
penetrate fairly far inside the halo, where they heat up to 
higher temperatures. The cold mode accretion ai z — 5.52, 
by contrast, is clearly directed along the intersecting fila- 
ments, and this coherent filamentary flow allows the galaxy 
to accrete gas from larger distances. A. Kravtsov (private 
communication) flnds similar behaviour in adaptive mesh 
simulations of an individual galaxy: the forming galaxy is 
surrounded by a hot gas halo, but cold filaments with co- 
herent flows p enetrate far inside this ha lo. The cluster scale 
simulations of iNaeai fc Kravtsov! ll2003fl do not include gas 
cooling, but they still show filamentary structures in which 
the gas entropy is far below that of the surrounding, spher- 
ically distributed gas. 

To quantify the evidence for filamentary cold accretion, 
we applied the following test. First, at the redshift output 



before the accretion output we find the radial vectors rf con- 
necting the centre of each resolved galaxy g to the positions 
of all particles i that it will accrete. Next we determine the 
normalised scalar product • r? = cos[rf,r^] for pairs of 
particles that accrete onto the same galaxy g. We compute 
the distribution of cos[rf ,r^] separately for all pairs of hot 
mode particles and all pairs of cold mode particles, summing 
the distributions for all galaxies g and normalising them to 
unity. Figure [TSl shows the results at z — 3 and z = 0. 

Hot mode particles show a nearly uniform distribution 
of cos[rf,r^], which is expected for a statistically isotropic 
accretion pattern. Cold mode particles, by contrast, show a 
strong peak at cos[r?, r?] « 1, indicating accretion from pre- 
ferred directions, i.e., along filamentary structures. The size 
of this peak underestimates the fraction of accretion that is 
filamentary, since particle pairs accreting from different fil- 
aments that feed the same galaxy produce cos[rf , r|] peaks 
at other angles, and averaging over the full galaxy popula- 
tion turns these multiple peaks into a uniform spread. Also, 
filaments are not perfectly straight, and they have a finite 
geometrical cross section, both of which tend to spread the 
cos[r^,r^] values. Nonetheless, this test shows a clear sta- 
tistical difference in the geometry of cold and hot accre- 
tion, consistent with the visual impression of Figure 1171 If 
we remake these plots separately for halos above and be- 
low Af = 2 X lO"Af0, we find that the directional signal 
in the cold accretion mode is stronger for high mass halos 
and weaker for low mass halos, while hot accretion remains 
approximately isotropic in both regimes. Some of the differ- 
ence between low and high mass halos could be numerical 
in origin, but we believe that it primarily reflects a more 
isotropic nature of cold accretion in low mass systems. 

From these clues, we can build the fo ll owing pic- 
ture of cold accre t ion. A s argued by iBinnevI (^73) and 
iBirnboim fc Dekell ij^QQ^)) accreting gas with short cooling 
times does not form a classic virial shock — its effective 
equation of state is too soft, and it passes through the virial 
radius ess entially in free fall. In th e spherically symmetric 
models of IBirnboim fc DekeJ (jioO^, there is a sharp tran- 
sition between low mass halos that do not have a virial 
shock and high mass halos that do, with the approximate 
condition being that a virial shock arises if the post-shock 
cooling time would be larger than the dynamical time. The 
important new ingredient highlighted by the simulations is 
the fllamentary nature of cold accretion. Departures from 
spherical symmetry allow galaxies to have a mix of cold 
an d hot accretion at one ti me, so the sharp transition of 
the lBirnboim fc DekeJ ll2003l) models is replaced by a steady 
change from cold mode domination at low mass to hot mode 
domination at high mass. A low mass galaxy has no virial 
shock anywhere, so it can accrete cold gas quasi-spherically, 
though fllaments still channel much of the flow. An interme- 
diate mass galaxy has a virial shock everywhere except at 
the points where dense filaments of cold gas penetrate the 
virial radius. These galaxies have a mix of cold, filamentary 
accretion and hot, quasi-spherical accretion. Gas in denser 
regions of the IGM is more likely to join a filament before 
accreting onto a galaxy, explaining the difference of density 
histories in Figure nm Finally, high mass galaxies have a 
virial shock everywhere and have only hot accretion. 

The characteristic halo mass at which we find equal 
cold and hot accretion rates, Af ~ 3 x 10^^ Mq, is a fac- 
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Figure 16. The average temperature (left panel) and density (right panel) evolution of cold mode (solid line) and hot mode (dashed 
lines) particles that are accreted onto galaxies aX z = 2. 



tor ~ 2 — 3 larger than t he transitio n mass calculated 
bv iBirnboim fc Dekell ^2003). However, 'Pcke l fc BirnboinJ 
i2005fl show that the derived transition mass depends on 
the precise criterion they use, in particular the radius at 
which they require a shock to form. For low metallicity gas, 
they find a transition mass ~ 1.5 x 10^^ for a shock at 
0.1-Rvir and ~ 10^^ A4q for a shock at Rvir, though the addi- 
tion of a UV background (as used here) would lower cooling 
rates and probably reduce these masses. Given the radical 
difference in calculational methods, the agreement between 
the SPH simulations and the spherically symmetric models 
is remarkably good, pointing to similar underlying physics. 

Semi-analytic models of galaxy formation distinguish 
between accretion with long post-shock cooling times, in 
which the cooling radius is smaller than the halo virial ra- 
dius, and accretion with rapid post-shock c ooling, in which 
the c ooling radius exceeds the virial radius ijWhite fc FrenkI 
I1991I) . In both cases, the cooling time is calculated assum- 
ing shock heating to the virial temperature, but in the latter 
case the gas is assumed to accrete onto the central galaxy 
at the free-fall rate instead of first forming a hot atmo- 
sphere in quasi-hydrostatic equilibrium. We have calculated 
the halo mass that marks the transition between these two 
cases for our adopted cosmology and cooling function, in- 
cluding the effects of the UV background. We calculate cool- 
ing rad ii assuming that halo gas follo ws an NFW density 
profile iNavarro. Frenk. fc Whitdfiggd) reduced by Q.b/ilm, 
and we take typica l concen trations as a function of redshift 
from lBullock et al.l (|200li) . If we define the cooling radius to 
be the radius at which the cooling time equals t{z), the age 
of the universe at redshift z, then we find a transition mass 
~ 7 X 10^° M0 at z > 1.0, roughly 2 to 4 times lower than the 
transition mass found in our simulations. At 2 < 1.0 there 
is no transition mass at all; in the presence of a UV back- 
ground, the cooling radius never exceeds the virial radius at 
the lower densities that prevail in low redshift halos. If we 
use either the freefall time or average dynamical time of the 
halo to define the cooling radius then the transition mass 
at high redshift is reduced by another factor of two to four, 
and there is no transition mass at z < 3 and z < 2.5. Thus, 
while the cooling radius transition in existing semi-analytic 



models bears a qualitative resemblance to the cold-hot tran- 
sition identified in this paper, a standard calculation of this 
transition yields a result that is quantitatively different and 
underestimates the significance of cold accretion (see the Ap- 
pendix for details of our calculation). The transition mass 
goes up if one ignores the UV background (thus raising the 
cooling rates) or assumes an gas density profile (thus 
raising the density at the virial radius) , but the first of these 
assumptions ignores a physical process present in the sim- 
ulation, and dark matter profiles, at least, are much better 
described by NFW profiles than by isothermal spheres. The 
cooling radius approach also fails to capture the filamentary 
nature of cold accretion, and because the cooling radius is 
calculated assuming virial temperature gas, it implicitly pre- 
dicts that most of the emission from even rapidly cooling gas 
will be at the X-ray wavelengths characteristic of these tem- 
peratures, while in the simulations this gas radiates much of 
its ener gy in atomic line t ransitions, hydrogen Lya in par- 
ticular iFardal et alJlioOlh . The exact amount of radiated 
Lya emission, however, needs further investigation, since it 
could be partially affected by the numerical broadening of 
shocks. 

One question remains; how does cold mode gas actu- 
ally enter a galaxy? Our ability to answer this question is 
limited because our simulations do not have good enough 
resolution at the relevant scales, but we know that gravita- 
tional infall to galaxy mass potential wells produces typical 
inflow velocities of ~ 100 — 300 km s~ ^ , and something must 
happen to the corresponding kinetic energy. (The simula- 
tion code conserves energy, so it cannot be lost to numerical 
effects.) We can identify three possibilities, all of which may 
operate to some degree. First, gas could penetrate inside 
the virial radius but be stopped by a strong shock close to 
the galaxy disk, heating to the virial temperature and con- 
verting its infall energy to thermal energy. The post-shock 
cooling time is then very short because the gas density near 
the disk is high. Some particles in our simulations exhibit ex- 
actly this behaviour, and A. Kravtsov's adaptive mesh sim- 
ulations also show shock heating far inside the virial radius 
at the terminations of cold filamentary flows (A. Kravtsov, 
private communication) . The lower left panel of Figure 1171 
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Figure 17. Accretion pattern around the largest galaxy in the L5. 5/128 simulation. Left hand panels show redshift z = 5.52, when the 
galaxy baryonic mass is 1.6 X IO^'^Mq and accretion is almost entirely dominated by cold mode. Green particles in each panel represent 
gas that will be accreted in cold mode by 2: = 5. Right hand panels show redshift z = 3.24, when the galaxy mass is 7.35 X IQ^^ Mq and 
accretion is almost entirely dominated by hot mode. Green particles in these panels show gas that will be accreted in hot mode by 2 = 3. 
Other particles in the top panels are colour coded by overdensity, from 0.67 (darkest blue) to > 2000 (white), in a region 4i?vir on a side 
(78/i~^ kpc and 202h~^ kpc in physical units at 2 = 5.52 and 3.24, respectively) and 2i?vir thick. Other particles in the middle panels 
are colour coded by temperature in the same region. Temperature scale starts from the lowest temperature (darkest blue) 1500 K on the 
left and 3000 K on the right panels, to the highest temperature (white); 5 X 10^ K on the left and lO'^ K on the right panels. Bottom 
panels have the same temperature colour-coding but zoom in to show a region IRvir on a side and 0.5iJvir thick. Vectors attached to 
particles show their projected velocities. 
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Figure 18. The distribution of tlie normalised scalar products of the radius vectors of accreting particles (see text) at the indicated 
redshifts. The dotted histogram is hot mode and the solid histogram is cold mode. 



shows some examples of such particles, with temperatures 
near Tvir in a quasi- spherical region of several kpc around 
the central galaxy. 

We formally count such particles as hot mode accretion, 
since they have a high Tmax, but physically they are more 
analogous to cold mode particles, since they do not shock 
heat near the virial radius. Because of our finite time resolu- 
tion, we miss the heating events for some of these particles 
and actually count them as cold mode. However, we think 
it is unlikely that most of the particles we identify as cold 
mode experience this kind of strong shock heating near the 
galaxy disk, in part because improving the time resolution 
makes only a small difference to the fraction of cold mode 
particles ( H4.2II a nd in part because of the cooling radiation 
arguments of Fa rdal et alJ 12001): the total amount of en- 
ergy radiated by accreting gas is of order the gravitational 
binding energy, and a large fraction of this energy emerges 
from low temperature gas. However, it is possible that high 
density accretion shocks near the forming disk are more com- 
mon in the real universe and are missed in our simulations 
because of their finite spatial resolution. 

A second possibility is that cold mode gas smoothly 
merges onto the galaxy disk like a stream of cars entering 
an expressway, converting its infall velocity to rotational ve- 
locity. This idea may seem outlandish at first, but the disk 
acquires its angular momentum from the accreting gas, so 
the required alignment of the disk with the accretion fiow 
may arise naturally. The large scale filamentary structure 
around galaxies persists for long periods of time, especially 
at high redshifts where the galaxies are highly biased with 
respect to the dark matter. A close examination of galaxies 
in our highest resolution simulation provides some evidence 
for this expressway mechanism, as one can see visually in 
the lower left panel of Figure 1171 but further analysis will 
be needed to assess its overall importance. In any event, the 
virial theorem suggests that this mechanism could provide 
only part of the energetic solution, since the galaxy's rota- 
tional energy should be smaller than its gravitational bind- 
ing energy. (Note, however, that the virial theorem cannot 
be applied with precision to non-isolated, time-dependent 
systems like forming galaxies, and the bookkeeping is fur- 



ther complicated by the fact that dark matter dominates the 
gravitational potential.) 

The third possibility is that infalling cold mode gas is 
slowed either adiabatically (i.e., by pressure gradients) or in 
a series of shocks that are individually too weak to cause 
strong heating. The left panel of Figure [T^ plots the average 
radial infall velocity of cold mode gas, scaled to the halo cir- 
cular velocity, as a function of r/vvir- We only consider cold 
mode gas that will be accreted onto the central galaxies by 
the next indicated redshift. We define the central galaxy as 
the most massive galaxy in the halo, which is almost always 
the galaxy closest to the halo center. We show results at 
z = 3 (dashed line, accreted by 2 = 2), z = 2 (solid line, ac- 
creted by z = 1) and z = 1 (dotted line, accreted by z = 0). 
Near the virial radius, the average infall velocity is similar to 
the halo circular velocity, but starting at r ~ 0.5-Rvir the gas 
decelerates smoothly, and by r ~ 0.1-Rvir the average infall 
velocity is only 10-20% of the circular velocity. The mean 
infall pattern is similar at all three redshifts, with the trend 
that the normalised velocities are lower at lower redshifts. 
There is a large scatter about the mean velocity, so some 
accreting particles must either shock heat near the disk or 
join the rotational flow as discussed above. However, based 
on Figure 1191 we tentatively conclude that deceleration by 
pressure gradients or weak shocks is the primary mechanism 
by which cold mode gas loses its infall velocity, at least in the 
simulations, radiating its dissipated energy as it slows. Ar- 
tificial cooling of gas in numerically broadened shocks (see 
M.'Sl could have some effect on this result, at least in its 
quantitative details. 

In contrast, the right panel of Figure ITUl shows the sit- 
uation for hot mode particles. We see that hot mode accre- 
tion more closely resembles the standard picture where the 
gas slows rapidly at a shock near the virial radius However, 
the zone where the gas slows down is significantly inside 
the virial radius at z = 3 (at about 0.6-Rvir) and moves 
steadily towards Rvii by z = 0. Furthermore, this zone is 
fairly broad, and the post-shock gas still has a mean ra- 
dial velocity that is a significant fraction of Vdrc, especially 
at high z. These complications refiect the departures from 
spherical symmetry and equilibrium discussed in i|6.1l Al- 
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Figure 19. Left: The median radial velocity of cold mode gas before it is accreted as a function of galactocentric radius. Radii are scaled 
to the halo virial radius, rvir, and velocities to the halo circular velocity, Hcirc- Results are shown at z = 1.0 (accreted by 2 = 0.0, dotted), 
2 = 2 (accreted by 2 = 1.0, solid), and 2 = 3 (accreted by 2 = 2.0, dashed). We only consider gas accreted onto the central galaxies of 
their parent dark matter halos. Right: Same for hot mode particles. 



though hot accretion shows no significant filamentary signal 
in Fig |18l from visual inspection we see that even in some of 
the biggest halos at high redshift a large fraction of the ac- 
creted gas comes from filaments. Filaments in these halos are 
hot (with temperatures well above the cold mode limit) but 
with temperatures much lower than Tvir- These filaments 
and mergers can deliver hot gas (as well as cold) well inside 
i?vir, with significant radial velocity. Halos at low redshift 
have better developed virial shocks and clearer separations 
from their filamentary surroundings, and the standard pic- 
ture becomes a better approximation. 

Further numerical investigations should significantly 
improve our phy s ical u nder standing of cold accretion. 
iBirnboim fc DekeJ ll2003h and lDekel fc BirnboimI J2005l) ar- 
gue that the physical criterion that determines whether a 
virial shock forms in their spherically symmetric calcula- 
tions is the ratio of the cooling time to the "compression 
time" under post-shock conditions. In future work, we will 
investigate whether this condition works as a local criterion 
for strong shock heating in our simulations. Our comparison 
of simulations with and without a UV background already 
points to an important role for cooling rates, and compar- 
ison of simulations with primordial and metal-line cooling 
and with difi'erent choices of fl;, and Q,m will help unravel the 
factors that determine the cold/hot transition mass. Most 
valuable will be very high resolution simulations that can ac- 
curately follow cold accretion streams all the way into galaxy 
disks, preferably performed with both SPH and adaptive 
mesh methods. 

6.3 Connection to Global Star Formation and 
Galaxy Structure 

As shown in Figure IT^ the cosmic star formation rate drops 
by a factor of ~ 11 from z = 3 to 2; = 0. Unfortunately, this 
star formation history cannot be directly compared to ob- 
servations because no one of our current simulations has the 
dynamic range to include all the galaxies that contribute sig- 
nificantly to the global star formation rate at all redshifts. 



In a forthcoming paper (Fardal et al., in preparation), we 
carefully combine the results from several simulations that 
span a range of volume sizes and resolutions to determine 
a global rate that can fairly be compared with the observa- 
tional results. We find that the global SFR is approximately 
constant from z = 6 to 2 = 2, then drops by a factor of 12 
from z = 2 to 2: = 0. The simulation predictions are con- 
sistent with the observations within their scatter, but this 
is not a stringent test because the observational estimates 
themselves span a wide range of values at any one redshift. 

Since the history of the cosmic star formation rate in 
our simulations closely tracks the history of gas accretion, as 
shown in Figure [T^ to understand the evolution of the cos- 
mic star formation rate it is sufficient to understand the evo- 
lution of the smooth gas accretion rate. From Figure|3]we see 
that both the cold and hot accretion rates decline at z < 2, 
with the cold mode decrease starting earlier and occurring 
more rapidly. As already discussed in JH] the transition in 
the global accretion rate from cold mode dominated at high 
redshift to hot mode dominated at low redshift simply re- 
fiects the increasing mass scale of galaxies and halos, since 
low mass systems are cold mode dominated and high mass 
systems are hot mode dominated at every redshift. However, 
the drivers of the declining accretion rate are somewhat dif- 
ferent in the two cases. 

First let us concentrate on the drop in the global cold 
mode accretion rate. In principle this could reflect a decrease 
in the space density of cold mode galaxies or a decrease in 
the amount of cold mode accretion per galaxy, or a combi- 
nation of the two. The rapid rise seen in Figure |4| reflects 
the increasing space densities of galaxies above our mass 
resolution threshold; most of these galaxies are cold mode 
dominated at high redshift. However, while the space den- 
sity of cold mode galaxies peaks at 2 ~ 2, it declines only 
very slowly thereafter. Most of the decline in cold mode ac- 
cretion must therefore be driven by the declining accretion 
rate per galaxy, as shown directly in Figure 1^ which plots 
the average gas accretion rate onto the central galaxy in ha- 
los of different masses as a function of redshift. For halos in 
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the mass range logMvir = 11 ± O.IM©, which is cold mode 
dominated at every redshift, the mean accretion rate per 
galaxy drops by a factor of ~ 12 between z = 2 and z = 0, 
accounting for the drop in the global cold mode accretion 
rate. 

The smooth solid curve in Figure I2UI shows the pre- 
diction of a simple analytical model along the lines of 
IWhite fc Frenk (1991): the accretion rate is just M = 
fgMvii where fg is the universal gas fraction (i.e., the uni- 
versal baryon fraction minus the fraction of baryonic mass 
in stars) and A/vir = (47r/3)-RvirPvir is the virial mass. The 
change in virial mass Mvir is determined through the cos- 
mological growth of the virial radius assuming an isothermal 
profile, fixed circular velocity and the cosmological evolution 
of pvir. For IO^^Mq halos, this model explains the simulated 
accretion rates quite well. One can view the steady decline 
of M as a consequence of the declining density at the virial 
radius, which is tied to the declining mean density of the 
Universe, or of the corresponding increase in characteristic 
dynamical times — the time for a halo to double its mass 
scales with the age of the Universe. 

At high redshifts the analytical model underpredicts the 
accretion rate by a factor ~ 2, for at least two reasons. First, 
much of the accretion is filamentary, allowing halos to draw 
material from beyond the virial radius and accrete more ef- 
ficiently than a spherical model predicts. Second, halos in 
this mass range are only ~ 2 times our resolution limit, 
and we are biased towards including halos that have higher 
accretion rates because we only count halos that have al- 
ready formed a resolved central galaxy. This effect operates 
at all redshifts, but it is more important at high redshifts 
when galaxies have had less time to grow, causing us to omit 
lower M galaxies from our average. In addition, our limited 
mass resolution may play some role in overestimating accre- 
tion rates. When we repeat the M analysis for the higher 
resolution, Lll/128 simulation at 2: = 3 and z = 4, the mean 
accretion rate is lower and agrees even better with the an- 
alytical model, though the rate is still higher by ~ 20% at 
z = 3 and ~ 40% at z = 4. 

Crosses connected by dotted and dashed lines in Fig- 
ure 1201 show the mean accretion rates for the central objects 
of more massive halos, with log Mvir /Af0 = 12 ± 0.15 and 
log Mvir /M0 = 13 ±0.15, respectively. Halos of lO^M© are 
dominated by hot mode accretion at all of the plotted red- 
shifts, and halos of 10^^ Mq are dominated by hot mode at 
all but the highest redshifts (see Figure |HJ. Unmarked lines 
show the corresponding predictions of the infall model de- 
scribed above; these correspond to the 10^^ Mq model curve 
(solid line) multiplied by factors of 10 and 100. In these 
higher mass halos, gas must cool before it can accrete onto 
the central galaxies, so it is not surprising that the pure in- 
fall model now overpredicts the galaxy accretion rates. How- 
ever, the infall model still captures the redshift dependence 
of these rates remarkably well, especially if one ignores the 
highest redshifts, where the number of 10^'^ A/0 halos is small 
and the W^^ Mq halos have a significant cold mode contri- 
bution. It is as if cooling just places a constant "tax" on 
the newly available halo gas, allowing an almost redshift in- 
dependent fraction of the gas to cool at a given halo mass. 
This fraction decreases from ~ 0.45 for IO^^A/q halos to 
~ 0.20 for W^^Mq halos, as one might expect given the 



1000 



" 100 



■S 10 





1013 


Mo 






- 


101= 

10" 


Mo 
M 

























_ -X ' 























0.2 0.4 0.6 

log(l+z) 



Figure 20. The average gas accretion rate onto the central 
galaxies of halos with logMvir/AfQ = 11 ± 0.1 (solid lines), 
logAfvir/A% = 12 ± 0.15 (dotted lines), and log Mvir /M© = 
13 ± 0.15 (dashed lines), as a function of redshift. Lines marked 
with crosses show results from the L22/128 simulation, while 
smooth, unmarked lines show the analytical infall models de- 
scribed in the text. 



higher virial temperatures and correspondingly longer cool- 
ing times in more massive halos. 

To model the cooling effects in simpler terms, we 
ha ve followed an analyt ical prescription along the lines 
of IWhite fc Freii^ il99ll) . assuming spherical halos with 
isothermal density and temperature profiles. We suppose 
that all gas that is able to cool accretes instantaneously 
onto the central galaxy. The accretion rate is, therefore, 
M = Anr'^^^ip^ooidr^ooi/dt, where pcooi is the gas density 
at cooling radius and the cooling radius, rcooi, is the radius 
where tcooi = tH and tn is the Hubble time. Therefore, the 
increase of the cooling radius in time defines the accretion 
rate. This model predicts accretion rates about 30%-60% 
lower than we find for lO'^^ and W^^ Mq halos, but it repro- 
duces the evolution of these rates fairly well. 

At least part of the difference in magnitude of accretion 
rates may be a result of ou r limited numerical resolution. 
ISpringel fc Hernauisd i2002l) find that the accretion rates 
can be overestimated in formulations of SPH like those that 
we use here when there are not enough gas particles in the 
galaxy halo. We have repeated the above analysis for the 
Lll/128 and L5. 5/128 simulations (at z = 3 and z = 4) and 
find that the accretion rates drop by 40%-50% for 10^^ Mq 
halos. However, each of these simulations has only one halo 
in this mass range, so this result is only suggestive. Fur- 
ther work using simulations with higher dynamic range will 
be needed to fully elucidate the mechanisms that drive the 
evolution of hot mode accretion. 

With these results in mind, we can build the following 
interpretation of the accretion rate evolution in Figure |^ 
and consequently of the evolution of the cosmic SFR in the 
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simulations. The accretion rate climbs very rapidly starting 
at 2: ~ 5 as more galaxies form above the resolution thresh- 
old — if we had higher resolution and a lower threshold, the 
accretion rate at these redshifts would be higher and would 
evolve more slowly. The number of galaxies in the low mass 
regime where cold accretion dominates reaches a maximum 
at z ~ 2 and declines slowly thereafter. The rapid drop in 
the global cold accretion rate, starting at 2 ~ 4, is driven 
mainly by the decreasing accretion rate per galaxy shown in 
Figure 1201 which reflects the increasing infall timescale in a 
lower density Universe. The mass threshold for hot accretion 
is higher, so the hot accretion starts to climb later and peaks 
at a lower redshift, 2 ~ 2. Furthermore, the number density 
of hot mode halos stays roughly constant at 2 < 2, and the 
average mass of these halos increases with time. (Note that 
higher mass halos always have higher mean accretion rates 
in both the simulation and the analytic model, even though 
the "cooling tax" makes the increase slower than M oc M.). 
The drop in global hot accretion is, therefore, also a conse- 
quence of the fast drop in the accretion rate per galaxy. This 
drop is slightly slower than for the cold mode owing to an 
increase in the characteristic mass of hot mode halos with 
time and to a slight increase in their number density. 

In broader terms, most accreting galaxies at high red- 
shift are low mass, do not have virial shocks, and connect 
directly to their large scale environment by filamentary "um- 
bilical cords." At intermediate redshifts, the galaxies expe- 
riencing the most accretion are intermediate mass systems 
with virial shocks penetrated by cool fllaments, accreting 
in a mix of hot and cold mode. At low redshift, there are 
many high mass galaxies with no cold accretion, and many 
low mass galaxies reside in high mass halos where they ex- 
perience primarily hot accretion. The filamentary structures 
that play such a vital role in cold mode accretion grow in 
scale as the Universe evolves, particularly around high den- 
sity regions. They eventually become too large to feed in- 
dividual galaxies and instead channel their flows into the 
hot IGM of groups and clusters. The combination of galaxy 
mass scales and large scale structure geometry drives the 
predicted dependence of the SFR on environment, allowing 
it to extend beyond the virial radii of large groups and clus- 
ters. 

One of the most striking features of the galaxy popula- 
tion to emerge from the SDSS is a rapid transition in typ- 
ical^ galaxy properties at a stellar mass M* ~ 3 x 10^" 
iKau ffmann ct al. 2003; Kannappan 2004). Galaxies below 
this mass tend to be actively star forming with lower sur- 
face densities, high gas fractions, and late type morpholo- 
gies, while galaxies above this mass tend to have old stel- 
lar populations, high surface density, low gas fraction, and 
early type morphology. The observed transition mass is re- 
markably close to the one at which we flnd a transition 
from cold mode domination to hot mode domination, so 
it is tempting to see the two as connected. It is not clear 
just how the second transition would bring about the first. 
However, if a galaxy accretes much of its mass along fila- 
mentary structures, this process is bound to influence its 
angular momentum, perhaps radically modifying the tra- 
ditional picture of disk formation in which the disk's spe- 
cific angular momentum distribution i s cl osely related to 
that of its parent dark matter halo iHo vlel 1194 91: iPeeblesI 
Il969l: iFaU and Efstathioul ll98ol: IMo. Mao, fc White! Il99a 



iBuUock etail l2001al : Ivan den BoschI 1200 iTl . In particular, 
cold accretion might be crucial to disk formation, with hot 
accretion contributing mainly to spheroid growth. The lower 
efficiency of hot accretion will also cause high mass galax- 
ies to have less star formation and older stellar populations, 
though it is not clear that this effect in itself will produce a 
transition as strong as the observed one. If, as we speculate 
in the next section, hot accretion is suppressed in the real 
universe by physical processes not represented in our simu- 
lations, then the transition would be much sharper. Further- 
more, galaxies substantially above the cold-to-hot transition 
mass would then have to be built by mergers, explaining 
their predominantly elliptical morphology. A tight connec- 
tion between cold accretion and disk formation would also 
help explain the observed morphology-density relation, since 
the cold accretion fraction depends strongly on environment. 
For the moment, these remarks are largely speculation, but 
we hope that future investigations of higher resolution simu- 
lations will deepen our understanding of the cold and hot ac- 
cretion processes and their connection to galaxy morpholo- 
gies and stellar populations. 

When the baryons within a dark matter halo cool to 
form the central galaxy in the standard hot mode picture, 
the dark matter densi ty near the galaxy incr eases through 
adiabatic contraction iBlumenthal et al]ll986ll . This results 
in an increase in the galaxy's circular velocity, which can 
lead to conflicts with observations. For example, it is difficult 
for models to simultaneously match the galaxy luminosity 
function and the zero-point of the TuUy-Fisher relation (e.g. 
Ivan den Bosch. Mo fc Yand boOSD . However, if the galaxy 
assembles through cold mode accretion, adiabatic contrac- 
tion might not occur or be greatly reduced since the dark 
matter and baryonic components assemble nearly simulta- 
neously. This might help alleviate some of the tension with 
observations, especially since most disk galaxies should be 
cold mode dominated given that on average they have lower 
masses (e.g. [Bell ct al. 2003). 

The most direct observational tests of cold accretion 
would come from detecting the cooling radiati on associated 
with g as infall onto galaxies. As emphasised bv lFardal et alJ 
i200ir) . the existence of an important cold accretion mode 
allows a significant fraction of the gravitational energy ac- 
quired by accreting gas to emerge in the Lya line, producing 
luminosities that are detectable with large telescopes. The 
main challenge to this test is separating the contribution of 
cooling radiation from the contribution of star formation, 
since galaxies with high accretion rates should also have 
high star formation rates. Because the Lya photons from 
cooling radiation are generally produced at larger galacto- 
centric radii, it should be possible to distinguish the two con- 
tributions, but accurate predictions of the angular and fre- 
quency distribution of Lya cooling radiation require radia- 
tive transfer calculations that have not yet been applied to 
3-d simulations. Filamentary gas flows that penetrate close 
to the galaxy disk and then shock should produce X-ray 
"hot spots" in the inner regions of some galaxy halos. Fur- 
ther investigation will be required to assess the observability 
of this phenomenon. High column densities of neutral hydro- 
gen near the disk might absorb much of the X-ray emission 
and re-radiate it in Lya (iBirnboim fc Dekelll2003f) . 
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6.4 Accretion and Feedback 

Although our simulations appear reasonably consistent with 
the observationally inferred history of cosmic star formation 
(Fardal et al., in preparation), they overpredict the z = 
baryonic mass function relative to the observational esti- 
mates of, e.g.,i2olc ct al. (2001) or Bell ot al. (200a) . Glob- 
ally, the predicted fraction of baryons in the form of stars 
and cold gas exceeds the observationally inferred fraction by 
a factor ~ 2 — 3. Of course, the integral of the mean cos- 
mic star formation rate must equal the mean stellar mass 
density (when corrected for the stellar mass loss and recy- 
cling of this gas), so the agreement with one and disagree- 
ment with the other implies that at least some of the ob- 
servational estimates of these quantities are themselves in 
conflict. However, the SFR estimates have substantial sys- 
tematic uncertainties, while the stellar mass densities appear 
reasonably secure if the assumed stellar IMFs are correct, so 
we think it more likely that the simulations are indeed pro- 
ducing overly massive galaxies. A related problem is that the 
simulations predict ongoing accretion and star formation in 
high mass galaxies, while in the real universe these systems 
are predominantly ellipticals exhibiting little or no recent 
star formation. Similar difficulties with the luminosities and 
colours of high mass galaxies appear in semi-analytic models 
of galaxy formation, unl ess they are specific ally modified to 
suppress them (see, e.g.. lBenson et al^l2003^ . 

Simply distinguishing between hot and cold accretion 
does not, of course, change the simulation predictions. How- 
ever, the conflicts mentioned above suggest that these sim- 
ulations are still missing some physical processes that play 
a significant role in galaxy formation, and these processes 
might have different effects on hot and cold accretion. One 
can imagine, for instance, that hot mode accretion might 
be more easily suppressed by some forms of feedback be- 
cause the hot mode gas is lower density and more isotrop- 
ically distributed. In the local universe, at least, it ap- 
pears that hot accretion is suppressed in groups and clus- 
ters, where X-ray emissivities imply short cooling times 
but X-ray spectra show that the predi cted gas at r<l 
keV is not present JPeterson et alJl2003ll . The leading hy- 
pothesis for explaining this conflict is that the central re- 
gions of the intracluster medium are heated by recurrent 
AGN activit y, thermal conduction from the outer regions, 
or both (e.g.,'Binnev & Tabor'l995';'Ciotti & Ostriker 200l|; 
[Naravan & Mcdvcdcv 2001; Ruszkowski & BoEolman 200^. 
The same processes might operate in lower mass halos at 
moderate redshifts. iMaller fc Bullockl ||2004) argue that a 
proper, multi-phase treatment of hot halo gas reduces pre- 
dicted accretion rates by a substantial factor, even with- 
out conduction or extra feedback, because the cooling gas 
forms dense clouds and the remaining gas has lower densi- 
ties and longer cooling times. This effect is omitted in stan- 
dard analytic treatments and would be missed by simula- 
tions (like ours) that do not resolve the ~ 10® A/0 scale of 
clouds formed by thermal instability. Whatever the mecha- 
nism, preferential suppression of hot mode accretion would 
lead to a sharper cutoff at the high end of the luminos- 
ity function, a more rapid decline in the cosmic SFR at 
2 < 1, and older stellar populations in massive galaxies. All 
of these changes would improve the agreement between our 
SPH simulations and observations, and they would also im- 



prove the agreement for most semi-analytic models. iBinnevI 
12004) has made similar points, drawing on our prel iminary 
results iKatz et alj| 2003|) and t he 1-d calculations of lBinnevI 
Jl977l) and lBirnboim fc Dekeil J2003fl f or theoretical motiva- 
tion (see also lDekel fc BirnboimlbOOstl . 

In our simulations, supernova feedback has little im- 
pact on galaxy masses because the supernova energy is usu- 
ally deposited in a dense medium, where it radiates away 
before it can drive a galactic wind. Simple models of a 
multi-phase interstellar medi um do not change this result 
JSpringel fc Hernauistll2003a^ . but collective effects may al- 
low real galaxies to drive winds more easily, and a number of 
groups have added ad hoc wind models to simulations (e.g., 
[Theuns et ai]l2002l : lAbadi et al.l20o3lSpringel fc Hernauisti 
l2003a[). The interplay of supernova-driven winds with accre- 
tion modes might be different from that of conduction or 
AGN feedback, because once a wind is launched the main 
challenge is getting it out of the galaxy halo, so that the gas 
does not simply return and form stars on a short timescale. 
In a galaxy fed by filamentary cold accretion, supernova- 
driven gas might encounter little resistance because there 
is no surrounding halo medium (see left panels of Fig I17II . 
In a galaxy dominated by hot accretion, on the other hand, 
the disk and bulge are surrounded by a shock-heated gas 
halo, and a wind may be trapped at small radius. Thus, one 
can easily imagine that expulsion of cooled gas by supernova 
feedback, in contrast to the suppression of cooling by AGN 
feedback or conduction, might operate more efficiently in 
low mass galaxies. This might provide a way of reducing the 
masses of low mass galaxies and flattening the faint end of 
the galaxy luminosity function, which would again improve 
agreement between the simulations and data. 

The conventional understanding of supernova feedback 
also suggests that it will be more effective in low mass sys- 
tems because gas is more easily driven out of shallow po- 
tential wells llDekel fc Silklll98d) . We have focused instead 
on the role of the surrounding hot gas, which we suspect 
is even more important, since simulations of the feedback 
process show that only low mass galaxies (Mgai^lO^M©) 
blow out significant amount s of gas if they are surro unded 
by conventional hot halos jFerrara fc Tolstovl 120001) . The 
form of the observed mass-metallicity relation can be under- 
stood if supernova-driven winds are efficien t below Mgai ~ 
SxIO^^Mq and inefficient at higher masses jTremonti et alJ 
I2OO4I) . This transition mass is close to the mass where cold 
accretion gives way to hot accretion in our simulated galax- 
ies, so it might emerge naturally if trapping by the hot gas 
halo is the correct explanation for it. 

The simulation predictions would also change in inter- 
esting ways if we assumed that cold accretion is associated 
with the convent ional stellar IMF (e.g.. Miller fc Scale 1973; 
iKennicutdllQSgTl but that hot accretion is associated with 
either a bottom-heavy IMF rich in brown dwarfs or a top- 
heavy IMF truncated below ^ 2Mq . In the first case, mas- 
sive galaxies would have high stellar mass-to-light ratios and 
therefore lower luminosities. In the second case, the mas- 
sive galaxies would have lower stellar masses (assuming that 
much of the recycled gas from evolved stars is driven from 
the galaxy), and they would have higher UV luminosities 
during their star-forming phases. Allowing IMF variations 
introduces a frightening level of freedom in galaxy forma- 
tion models, but with two physically distinct accretion pro- 
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cesses one should at least consider the possibility that the 
associated star formation is different. 

For the moment, these proposals are largely speculation. 
We are carrying out further investigations to see how the 
predicted properties of the galaxy population would change 
if we assume preferential suppression of hot mode accretion 
or different stellar IMFs for the two modes. Such adjust- 
ments may make distinctive observational predictions that 
could compensate for the extra freedom they introduce, and 
their physical underpinnings can be investigated with higher 
resolution simulations focused on these questions. Given the 
variety of possible assumptions, semi-analytic calculations 
might be a sensible approach to exploring models that treat 
cold and hot accretion differently. 



7 CONCLUSIONS 

In SPH simulations of galaxy formation, a substantial frac- 
tion of the gas accreted by galaxies is never heated close to 
the virial temperature of the galaxy potential wells. The im- 
portance of this cold accretion mode, relative to the hot ac- 
cretion mode envisioned in the traditional picture of galaxy 
formation, depends strongly on galaxy mass. Most galaxies 
below a baryonic mass Mgai ~ 10^" '^ Mq or dark halo mass 
Afhaio ~ lO^^ ^M© accrete primarily in cold mode, while 
more massive galaxies accrete primarily in hot mode. Since 
high mass galaxies are built from lower mass systems, even 
galaxies at the top end of the luminosity function today ac- 
quired a significant fraction of their mass via cold accretion. 
Globally, the increasing mass scale of galaxies means that 
cold accretion dominates at high redshift and hot accretion 
at low redshift. The ratio of cold to hot accretion is also en- 
vironment dependent, mostly because higher mass galaxies 
are more common in dense environments, and partly be- 
cause low mass systems in dense environments have a larger 
hot accretion fraction than their isolated counterparts. Hot 
mode accretion is quasi-spherical, but cold mode accretion 
is often directed along filamentary channels, increasing its 
efficiency. 

A number of arguments suggest that these are phys- 
ical results, not numerical artifacts. A quantitatively im- 
portant cold mode appears in simulations covering a fac- 
tor of 512 in mass resolution, though the lowest resolu- 
tion simulations do not include a photoionizing UV back- 
ground and therefore predict a noticeably different temper- 
ature distribution for cold mode gas. The transition mass of 
Afgai ~ lO^^'^M© is stable in simulations covering a factor 
of 64 in mass resolution. Cold accretion is seen in other SPIf 
simulations, including one with a substantially different im- 
plementation of shock heating. Cold filamentary flows are 
found in adaptive mesh simulations, though these cannot 
yet test whether gas gets into galaxies without experiencing 
strong shock heating. Finally, the transition mass in our sim- 
ulations agrees (to within a factor of ~ 2—3) with that found 
by^rnboim & Dokcl (2003) using entirely independent an- 
alytic methods and high-resolution 1-d calculations. Despite 
this litany of tests, we cannot rule out the possibility that 
cold accretion is an artifact of finite resolution of shocks in 
SPH simulations. Confirmation with still higher resolution 
simulations and with hydrodynamic methods that directly 
impose shock jump conditions at velocity discontinuities is 



therefore desirable. The transition mass separating hot and 
cold accretion is quantitatively different from the transition 
mass separating rapid and slow post-shock cooling in semi- 
analytic models of galaxy formation, at least as calculated 
by our implementation, which we believe to be the standard 
one (see H6.2I and the Appendix). However, semi-analytic 
models could be easily revised to incorporate a cold mode of 
gas accretion calibrated on hydrodynamic simulations like 
those analyzed here. 

Cold accretion could have important implications for 
the cosmic star formation history and for the correlations 
of galaxy star formation rates, stellar populations, and mor- 
phology with environment. If cold and hot accretion are af- 
fected differently by stellar or AGN feedback, or if they are 
associated with different stellar IMFs, then the simulation 
predictions would change in interesting ways, perhaps re- 
solving conflicts with the observed colours of ellipticals and 
with the observationally inferred galaxy baryonic mass func- 
tion. Our predicted transition mass is close to the observed 
transition mass at which ty pical galaxy properties change 
markedly iKauffmann et a'Ll l2003). and we have speculated 
on possible connections between the two. In particular, if 
AGN feedback, multi-phase cloud production, or some other 
mechanism strongly suppresses hot accretion in the real uni- 
verse, then galaxies above the transition mass would have 
little recent accretion and star formation, and they would 
have to be built by mergers of lower mass systems, which 
would produce spheroidal morphologies. Many questions re- 
main about both the underlying physics and the observa- 
tional implications of cold accretion, but we hope that higher 
resolution simulations and comparisons to data will answer 
these questions over the next few years and solidify our un- 
derstanding of how galaxies get their gas. 
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APPENDIX A: THE RAPID COOLING 
TRANSITION IN SEMI-ANALYTIC 
CALCULATIONS 

Semi-analytic models of galaxy formation, which trace back 
to the formalism introduced bv lWhite fc Fren^ il99 J) , cal- 
culate the growth of galaxies fed by the accretion of gas 
cooling within dark matter halos. While cooling rates are 
calculated assuming that the gas initially heats to the halo 
virial temperature, these calculations do distinguish between 
a regime of slow cooling, in which accretion from a quasi- 
hydrostatic gas halo is regulated by the cooling rate, and a 
regime of rapid cooling, in which no quasi-hydrostatic at- 
mosphere forms and gas accretion is regulated by the infall 
rate. Given the qualit ative resemblance of this criterion to 
the one introduced bv lBirnboim fc Dekell ijioO^, it is inter- 
esting to ask whether the predicted transition between infall 
and cooling dominated accretion is quantitatively similar to 
the transition between cold and hot accretion found in our 
simulations. In this Appendix, we describe our calculation of 
the infall-cooling transition mass, discussed earlier in H6.2I 
We assume that galaxy halos initially contain the uni- 
versal baryonic fraction of gas, Q.i,/Q,m, with a constant gas 
temperature equal to Tvir. We m odel the initial gas density 
distribution as an NFW profile iNavarro. Frenk. fc Whita 
Il996l) . and, for simplicity, assume a constant concentration 
for the NFW profile s, which is typical for 10^^ halos 
f Bu llock et ai]l2001hD : c = i^vir/r^ = 15/(1 + z). i?vir is the 
virial radius and is the scale radius of the NFW halo. 
For a given halo virial mass, Mvir, -Rvir is defined as the ra- 
dius within which the enclosed density equals Avir times the 
mean mass density of the universe at that redshift, i.e. 

1/3 



-Rvir — 



Here 



/ 3 Mvir 
\^47r Avir/Om 



(Al) 



the 



iBrvan fc NormanI lll998D. 

ioj lie 



approximation for Avir from 
which gives similar results 
199a). For comparison we also 



to iKita^Mnafc Suto 
consider an alternative model where the gas density follows 
that of a singu lar isothermal sphere ( SIS) with Tgas = Tvir. 

FoUowing IWhite fc FrenkI diggiTl . we define the local 
cooling time of the halo gas, tcooi, as 

3 ksTpgir) 



tcooi{r) 



(A2) 



2 fg^J.mJ,nl{r)A{T) ' 

where fcs is the Boltzmann constant, fi is mean "molecu- 
lar" weight, which we assume to be that of a fully ionized 
gas, fg is the fraction of halo gas available for cooling, irip 
is the proton mass, uh is the hydrogen density, pg is the 
total gas density, and A(T) is the cooling function. Like in 
the simulation, we assume that the gas is primordial with a 
hydrogen mass fraction of 0.76. The cooling time determines 
the cooling radius rcooi at which gas can radiate its therma l 
energy in a specified amount of time. lwhite fc FrenkI lll99j) 
define the cooling radius by setting fcooi(?'cooi) equal to the 
Hubble time t{z), the age of the Universe at redshift z. How- 
ever, difi'erent groups use different characteristic times in 
their semi-analytic calculations. Therefore, to cover a range 
of possibilities we compare the cooling time to three time 
scales: the Hubble time t{z), the dynamical time defined as 



^dyn 



3tv 
16Gp{z) 



(A3) 
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where p is the average enclosed density, and the free fall 
time, which is = tdyn/V^- Therefore, the cooling radius 
is define d as the radius where t cooi(rcooi) = t{z), tdyn(z), or 
iff (2). In lWhite fc Frenld il99ll) . infall dominated accretion 
takes place when rcooi > Rvii, and cooling dominated accre- 
tion takes place when rcooi > Rvir- The transition halo mass 
between the two regimes is thus the halo mass for which 

^cool — Rvir • 

To fairly compare this mass with the transition mass 
between cold and hot mode accretion in our simulations, 
we use the same cosmology (see § 12. H and the same cool- 
ing function (see KWH96) as in our simulations, neglecting 
only inverse Compton cooling. Since some gas has already 
transformed into stars, the fraction of the gas available for 
cooling, fg, is less than one and decreases with redshift. We 
approximate fg = i~ 0.28/(1 + z), which roughly matches 
the global fraction of the baryons in the simulation that are 
not in stars or cold galactic gas at a given redshift. 

The UV background affects gas cooling in our simula- 
tions after z = 6, so to make a fair comparison it is necessary 
to include it in our calculation. As mentioned in § 14.41 for 
typical densities at the virial radius, the UV background 
can significantly alter the hydrogen and helium line cooling 
peaks in the cooling function, making the cooling times at 
the virial radius siginficantly longer. 

In Figure lAll solid lines show the transition mass be- 
tween infall and cooling dominated accretion for an NEW 
density profile in the presence of the UV background. 
Squares correspond to a definition of rcooi with t^ooi ~ i(-z), 
triangles to tcooi ~ tdyn{z), and the star to tcooi ~ tff(z). 
For direct comparison we plot the mass where the transi- 
tion between cold and hot mode occurs in our simulations 
as the dot-dashed line with circles. At z > 1 the transition 
mass in the simulation is a factor of 2 — 4 higher than the 
highest of the semi-analytic transition masses, defined for 
tcooi{Rvii) ~ tiz)- At z < 1, there is no halo mass for which 
tcooi(^vir) < i(-z), so the semi-analytic calculation predicts 
no infall dominated halos at all. If we define the transition 
with respect to tdyn or tg, the semi-analytic transition mass 
drops by another factor of 3 — 4 at high redshift, and infall 
dominated halos disappear after z — 2.5 (tdyn) or z — 3 (tff). 

The long-dashed curve in Figure IXTl shows the efi'ect of 
omitting the UV background, which increases cooling rates 
and therefore increases the transition mass. With no UV 
background, there are infall dominated halos all the way 
to 2 = 0, though the transition mass is still lower than 
the cold/hot transition mass in our simulations with a UV 
background. Furthermore, as discussed in H4.4I the increased 
cooling in the absence of a UV background leads to a sub- 
stantially higher transition mass in our simulations them- 
selves. Since our simulations without a UV background are 
lower resolution or run only to 2: = 3, we do not have a 
very accurate measurement of this increase, but it appears 
to be at least a factor of four. The dotted curve in Figure IXTl 
shows the effect of assuming an SIS gas profile instead of an 
NFW profile, again with no UV background. Since an SIS 
profile has a higher gas density at the virial radius than an 
NFW profile, the transition mass increases. However, while 
gas and dark matter profiles inside halos will generally be 
different to some degree, the NFW assumption seems likely 
to be more realistic ne ar the virial radius. For co mparison we 
also plot the model of lBirnboim fc DekeJ ll2003ll (zero metal- 



licity line from the upper panel of their Figure 11), which 
is motivated by their ID simulations. Their cosmology is 
different from ours, but the baryonic mass fraction, one of 
the main fa ctors determining the transition mass, is similar. 
iBirnboim fc Dckol (2003) do not include a UV background 
in their calculations, and their model agrees fairly well with 
the "NFW without UV" fine in Figure FCT though the red- 
shift d ependence is somewhat different. IPekel fc BirnboimI 
i200^ discuss the dependence of the derived transition mass 
on details of the calculational assumptions, gas metallicity, 
and cosmological parameters. 

We conclude that the conventional definition of the 
transition between infall and cooling dominated accretion 
does not correspond quantitatively to the transition between 
cold and hot mode accretion found in our simulations, if one 
makes the same physical assumptions in the semi-analytic 
calculation. However, it may be possible to alter the semi- 
analytic calculation, for example by defining the transition 
at some fraction of the virial radius instead of at i?vir itself, 
so that the correspondence is better. It would be partic- 
ularly interesting to see whether a definition of transition 
mass calibrated to match hydrodynamic simulations for one 
set of cosmological parameters and cooling rates then repro- 
duces the results of simulations that assume different pa- 
rameters or different cooling rates (e.g., because of different 
metallicity or a change of the UV background) . A successful 
analytic model of this sort would provide physical insight 
into the distinction between cold and hot accretion, and it 
would be a useful practical tool for investigating the effects 
of parameter variations. 
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Figure Al. Transition halo mass between infall and cooling 
dominated accretion, calculated in a manner similar to that of 
IWhite fc FrenlJ I1991I) . Squares, triangles, and the star at 2 = 3 
assume an NFW gas density profile with a UV baclsground and 
a cooling radius definition tcool ('"cool) = ^(-2)1 *dyn(2), and tg{z), 
respectively. With a UV background, all halos are cooling dom- 
inated below 2 = 1, 2.5, or 3, respectively, for these three cases. 
The dashed line assumes an NFW profile but no UV background, 
and the dotted line assumes an SIS profile with no UV back- 
ground, both for tcool(''cool) = ^(z). Filled circles show the tran- 
sition halo mass at which we find equal amounts of cold and hot 
accretion in the L22/128 simulation, which includes a UV back- 
ground. For comparison we also show results from the model of 
iBirnboim fc DekeJ i20n3h (see text). 
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